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ABSTRACT 

Accurate fixed-frequency virtual height data collected on a digital 
ionosonde is analyzed to measure speed and direction of traveling iono- 
spheric disturbances by matching the experimental data with virtual height 
and echo amplitude obtained from a simple model of the disturbed ionosphere. 
Several data records analyzed in this manner indicate speeds of 400 
to 680 msec -1 and a direction of propagation from north to south. The 
digital ionosonde collects virtual height data with a time resolution 
of 10 sec and a height resolution of less than 300 m. 
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1 . INTRODUCTION 


Virtual height recorders were among the first instruments used 
to study traveling ionospheric disturbances. Virtual height data are 
generally recorded on photographic film with a height resolution on the 
order of a few km. With the availability of high speed computers for 
real time data collection, better resolution can be obtained by digitizing 
the output of a virtual height recorder and analyzing the data with 
a computer. Besides better height resolution, computer algorithms provide 
noise reduction and greater flexibility in collecting and recording the 
data. 

1.1 Statement of the Problem 

The purpose of this work is to study traveling ionospheric 
disturbances (TIDs) i n the F region by comparing accurate fixed- 
frequency virtual -height data from a digital ionosonde with 
virtual heights calculated from simple models of the disturbances. 

Horizontal speed, wavelength and direction of the disturbances 
are calculated from virtual -height data. The ionosphere is modeled 
as a curved reflector of radio waves. Virtual height and focusing 
effects for a reflector are calculated and compared with experimental 
data. The shape and size of the reflector are varied until the 
calculated height and amplitude match the experimental data. The 
horizontal wavelength of the disturbance is the horizontal size of 
the reflector and the speed of the disturbance is the wavelength 
divided by the period of the disturbance. 

To determine the direction in which the disturbance is 
moving, the separation between the ordinary and extraordinary 



wave reflection points is calculated. The time delay between 
the appearance of prominent features of the disturbance at the 
two reflection points is measured to determine the component of 
speed projected in the N-S direction. The angle at which the 
disturbance is propagating relative to North is the arc -cosine of the 
N-S projection of the velocity divided by the speed of the disturbance. 
The direction measured is either east or west of North. The 
ambiguity cannot be resolved with the present equipment. 

1.2 Outline of the Thesis 

Chapter 2 reviews gravity wave and radio wave propagation 
theory and describes how radio waves reflected from the ionosphere 
are used to study gravity waves. 

Chapter 3 describes the data collection system and contains 
an analysis of the accuracy of the algorithm used to determine 
virtual height. 

Chapter 4 discusses the results of the analysis of data 
by methods described in Chapter 2. 

Chapter 5 contains the conclusions and suggestions for further 

work. 

Operating procedures and listings of the programs for collecting 
and analyzing virtual height data are given in the appendices. 
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2. ANALYSIS OF TRAVELING IONOSPHERIC DISTURBANCES BY 
VIRTUAL- HEIGHT MEASUREMENTS 

Traveling disturbances in the F region have been observed for more than 
20 years. Munro [1950] used fixed- frequency virtual -height measurements to 
determine speed, direction and extent of these disturbances. Bines [1959] 
developed a theory of gravity waves to explain the disturbances observed by 
Munro and others. 

In this chapter the effects of gravity waves on the propagation of 
radio waves are discussed and also how radio waves are used to measure the 
properties of traveling ionospheric distrubances . 

In Section 2.1 the properties of gravity waves and their effects on 
the ionosphere are discussed. 

In Section 2.2 the propagation of radio waves which retrace their path 
after reflection from a disturbed ionosphere is discussed. 

In the following sections, techniques are developed to measure the 
speed and direction of the disturbance from properties of radio propagation 
in an ionosphere disturbed by gravity waves. 

2 . 1 Gravity Waves 

Sound waves propagating through the ionosphere are called gravity waves 
because of the effect gravity has on them. Gravity causes the density of the 
atmosphere to decrease with height. This horizontal stratification causes 
the propagation to be anisotropic. 

A derivation of the wave equations can be found in Yeh and Liu [1972]. 
The dispersion relation for a gravity wave in an isothermal atmosphere is 
given by 
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where n - k fk 
x x o 

n = k /k 
z s o 


c s ~ speed of sound 

W a = Y 5 r / 2 < 2 g acoustic cutoff frequency 


1/2 

u. = (y-1) g/o buoyant or Brunt frequency. 

u S 

The coordinates are set up with z vertical, and with k Q in the x-s 

plane 

^ = \T<- q \ = tl 16 wave number 

k and k are the components of k in the x and z directions 
x z o 

respectively. 

y = the ratio of specific heats 
g - the gravitational acceleration 
For frequencies less than the Brunt frequency, the refractive index 
curve is shown in Figure 2.1. For these waves, called internal gravity 
waves, the phase velocity is directed downward when the waves are traveling 
upward. 

Waves with frequencies greater than the acoustic cutoff frequency are 
called acoustic gravity waves and have an elliptical refractive index curve 
as shown in Figure 2.2. 

For frequencies between the Brunt and acoustic cutoff frequencies, 
either k or k is imaginary. Freely propagating waves cannot have fre- 

CC 

quencies in this range, called the cutoff region. 




INTERNAL GRAVITY WAVES 


Figure 2.1 Refractive index curve for an internal gravity 
wave from Georges [1967] . 





ACOUSTIC WAVES 

Figure 2,2 Refractive index curve for an acoustic gravity 
wave from Georges [1967] . 
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Typical values for co^ and to^ at 250 km are .483 and .436 rad sec 
corresponding to periods of 13 and 14.4 min respectively. 

Associated with the propagation of a gravity wave is the motion of 
small air parcels which move in elliptical orbits. Above the dynamo region, 
a force applied on the ions by collisions with neutral particles tends to 
drag them along with the air parcel but because of the earth's magnetic 
field, the ions are constrained to move in a spiral path along the field 
lines. Neglecting diffusion, the resulting ion velocity is equal to the 
component of the air parcel velocity parallel to the magnetic field. 
Electrons move with the ions to maintain neutrality. 

Figure 2.3 shows the iso-ionic contours for a traveling ionospheric 
disturbance observed at Arecibo by Harper [1972], In the next section, the 
propagation of radio waves in a disturbed ionosphere such as shown in 
Figure 2.3 is discussed. 

2.2 Radio Waves as a Tool for Studying Traveling Ionospheric Disturbances 

Almost all data on traveling ionospheric disturbances comes from radio 
propagation experiments. The earliest techniques relied on the time delay 
of waves reflected from the ionosphere. A description of wave propagation 
in a horizontally stratified ionosphere is given in Budden [1961]. 

For radio waves in the F region with frequencies greater than 1 MHz, 
the refractive index is given by 


v = 1 


1/2 


-t t (r 2 +y L 2 ) 1/2 


( 2 . 2 ) 


The axis are set up with z vertical and the magnetic intensity vector B 
in the x-z plane. Z and n are the x and z direction cosines of 2. 
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Figure 2.3 Iso-ionic contours in the presence of a gravity wave 
from Harper [1972] . 
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T = Y t 2 /2(1-X) 

y l = nY 
? T = W 

Y = = t S e/m^ | gyrofrequency 

2 2 2 

X = (<*> n A0 £n)^ = He ItJHg plasma frequency 

is the permittivity of free space 
N is the electron density 
e is the charge on an electron 

m is the mass of an electron 

e 

The refractive index has two values corresponding to the ordinary 
and extraordinary polarizations. For a horizontally stratified ionosphere, 
a wave transmitted vertically and reflected from the ionosphere will re- 
trace its path and return to the ground at the point where it originated. 

Reflection occurs where y = 0. The ordinary wave path is deflected towards 

2 2 

the magnetic pole and reflected when X = 1 or N = /® • The extra- 

ordinary wave path is deflected towards the equator and reflected when 

x = i - y. 

Gravity waves cause the ionosphere to be nonhorizontal ly stratified 
and waves which retrace their path are not vertically incident on the iono- 
sphere. Under these conditions, the reflection conditions X = 1 and X = 

1 - Y for the two polarizations still hold, but the path is deflected in 
the x and y directions from the path it would take in a horizontally strat- 
ified ionosphere. 

The most convenient quantity to measure, relating to disturbances in 
the ionosphere, is virtual height h 1 , the equivalent height of the reflec- 
tion point if the wave traveled at the speed of light a. If the time delay 
between the transmission of a wave and its reception at the same point after 
being reflected from the ionosphere is T, the virtual height is h * = oT/2. 
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For the purpose of studying gravity waves, the ionosphere is modeled 
as a curved reflector with radio waves that retrace their path reflecting 
at right angles to the surface as shown in Figure 2.4. The reflector 
moves horizontally with the same velocity as the gravity wave. The shape is 
assumed to be constant with time and not to vary in the y direction. The 
virtual height of a wave reflected from point A is the slant height from the 
ground to point B on the reflector. The virtual height measured from a 
fixed location will vary as the reflector moves horizontally with velocity v. 
To calculate virtual height for a reflector, it is easier to think of the 
reflector remaining stationary and the transmitter and receiver moving along 
the ground with velocity -v. 

Virtual height measurements from the digital ionosonde provide the 
period and a measure of the amplitude of the disturbance directly. Several 
methods for determining horizontal wavelength, speed and direction are dis- 
cussed in the following sections. 

2.3 Focusing of Radio Waves due to Ionospheric Irregularities 

One property of a curved reflector is to focus or defocus the waves re- 
flected from it. The intensity of the radio waves reflected back to the 
receiver is related to the curvature of the reflecting surface. Whitehead 
[1956] derives a formula to relate the amplitude of the received wave to 
the curvature of the reflector using the approximation that the transmitter 
illuminates the reflector with constant power per unit area. The following 
derivation is similar but drops that assumption. Figure 2.5 shows the geo- 
metric relation between the reflector and the transmitter. At the point of 
reflection, the reflector has a one dimensional radius of curvature R. The 
amount of power illuminating the reflector over an area of length J and 
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Figure 2. 



5 Focusing of radio waves by a curved reflector. 
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unit width is focused onto an area of length L and two units wide. The 
amount of power illuminating the reflector per unit area is proportional to 
h . The mirror formula is 1 /h + 1/s = 2/H. By similar triangles 
J/L = s/{s-h ) = R/2(R-K) . The amplitude squared of the returning wave is 
proportional to the power per unit area at the ground which is equal to the 
power per unit area at the reflector times the focusing factor R/2(R-h). 


A 


2 



R 

R-h 


(2.3) 


If the shape of the reflector is given by H(x) , the radius of curvature 
R is given by 


. i * iMzL 2 3/2 iiS^L (2.4) 

A dx 2 

Focusing depends on the horizontal structure of the reflector so the 
amount of focusing can be used to determine the reflector shape. This is 
discussed further in the next section. 

2.4 Effects of More Than One Reflection Point 

When the center of curvature of the reflector is above the ground, 
three rays reflected from the dome-shaped reflector retrace their paths 
bach to their source simultaneously as shown in Figure 2.6, 

For an assumed reflector shape H(x), the virtual height of a ray re- 
flected from a point x q on the reflector is given by 

h ' = H(x q ) (1 +m 2 ) 1/2 


(2.5) 
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where m is the slope of the reflector at the reflection point and the 
point on the ground where the ray originated is given by 

x - x q + H[x o )m (2.6) 

A program for a Hewlett-Packard Model 9830 calculator plots virtual 
height and amplitude variations caused by focusing for any reflector shape. 
Operation of the program is discussed in Appendix III. The average height, 
peak-to-peak amplitude, horizontal size and shape of the reflector are ad- 
justable. The reflector shape, virtual height and amplitude are plotted 
on a Hewlett-Packard Model 9862 plotter. An example is shown in Figure 2.7. 
The horizontal axis is in km. The graph can be thought of as the height 
and amplitude for a stationary reflector as the transmitter and receiver are 
moved in the horizontal direction. If the reflector is moving and the trans- 
mitter and receiver are stationary, the horizontal axis is proportional to 
time. The proportionality factor is the speed of the reflector. 

If the virtual height and amplitude calculated from an assumed re- 
flector shape matches the experimental data, the horizontal length of the 
gravity wave is equal to the length of the reflector. The horizontal phase 
speed of the disturbance is the wavelength divided by the period of the 
disturbance. 

2.5 Second Reflections 

A ray which travels from the ground to the ionosphere and back twice 
is called a second reflection or a two-hop ray. The simplest two-hop ray 
path is shown in Figure 2.8. The two-hop ray is not incident at right 
angles on the reflector and does not retrace its path after being reflected 



lIGHT (km) 





Figure 2.8 One- and two-hop ray paths 
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from the ionosphere as the one-hop ray does. Because of this, the two-hop 
ray is reflected below the height where X = 1 and X = 1 - X for the ordinary 
and extraordinary polarizations respectively. For the reflector slopes of 
less than 20 degrees normally found in the ionosphere, the two -hop ray is 
nearly incident at right angles to the reflector and the one-hop reflector 
is used as an approximation of the reflection height of the two-hop ray. 

The virtual height and amplitude for a two-hop ray can be calculated 
for a known reflector tf(a :) in the same manner as for the one-hop ray. The 
virtual height for a ray which is reflected at x q is 

h' - H(x q ) (1+1/cos (2 tan' 1 (m))) (2.7) 

where m is the slope of the reflector at x q . The ray originates at 

x = x q + tan (2 tan _1 (m)) (2.8) 

Focusing is calculated from the geometry shown in Figure 2.9. For 
simplicity, instead of being reflected at the ground, the ray passes through 
and is reflected from a mirror ionosphere on the other side of the ground 
line. From the mirror equation l/e + l/w~ 2/R and 1/t + 1/(2 H^x^-w) = 2/R. 
By similar triangles K/J - (2 H(jx^~u))fw and L/R = (t-s)/t. The amplitude of 
the ray as recieved is equal to the power per unit area along J times the 
focusing factor J/L. 


A 2 = L. —JL ±- 

2 t-s 


(2.9) 




Figure 2.9 Geometry of two-hop focusing. 
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The program which calculates the height and amplitude for a one-hop 
ray path, also calculates the two-hop virtual height and amplitude. Figure 
2.10 shows the two-hop virtual height and amplitude for the reflector shape 
used in Figure 2.7. The two-hop ray path data are used in the same manner 
as the one-hop data to determine the shape and size of the reflector. 

A method of determining the direction of the disturbances is dis- 
cussed in the next section. 

2.6 Lateral Deviation 

For a horizontally stratified ionosphere, the deviation of the wave 
path from vertical, the lateral deviation, is given by Budden [1961] as 


dx _ l n y 2 (1-iQ 2 

“ C4j l 2 (i-z) 2 + y/) 1/2 


( 2 . 10 ) 


The - sign applies to the ordinary polarization and the + sign applies 
to the extraordinary polarization. The equation is integrated numerically 
to give lateral deviation. Table 2.1 gives the lateral deviation versus 
height above the bottom of the ionosphere for the two polarizations. For 

this calculation the F region of the ionosphere is modeled as a parabola 

. 12 -3 
with a width of 400 km and a maximum electron density of 10 m 

For an ionosphere disturbed by a gravity wave, the reflection points 
oscillate around the points for an undisturbed ionosphere and the separa- 
tions calculated are average separations. 

Since the reflection points are separated in the horizontal and verti- 
cal direction, the time delay between the appearance of prominent features 
of the disturbance at the two reflection points depends on the vertical 
phase variation and the horizontal velocity of the disturbance. The phase 



HEIGHT (km 



HORIZONTAL DISTANCE (km) 

Figure 2.10 Virtual height and amplitude for the two-hop ray calculated for the 
reflector in Figure 2.7. 
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progression for disturbances with periods greater than 25 min is downward. 
Georges [1967] measured the angle of the phase front for a disturbance with 
a period of 40 min at 83° from vertical. For shorter periods the angle of 
the phase front decreases to 45°. 

Figures 2.11a, b, and c indicate the reflection points and the phase 
front of the distrubance looking west, down and perpendicular to the velocity 
vector of the disturbance respectively. If the maximum for the ordinary 
reflector is at the ordinary reflection point, the maximum for the extra- 
ordinary reflector has a distance of D + D tanif) to go to reach the extra- 
ordinary reflection point. If the time it takes foT this to happen is t* 
the projection of the velocity in the southerly direction is 

V = (jD + D tan<j> secQj/x. The cosine of the angle 0 is the speed of the 
s cc a 

reflector v divided by the velocity projection in the southerly direction. 

cose = V/V = Vx/(D +D tan<J> sec0) (2. 11) 


This can be solved for 0 

Vt - D tan<j> 

0 = cos" 1 f (2.12) 

x 

The angle 0 is measured either east or west of north. The ambiguity cannot 
be resolved with the present equipment. 

The velocity v is determined by methods described in the last three 
sections. Chapter 4 contains the results of applying the methods of this 
chapter to data collected by the digital ionosonde described in Chapter 3. 





(a) LOOKING WEST 


(b) LOOKING DOWN 


D z cot<£ 



(c)LOOKING NORMAL TO DIRECTION 
OF GRAVITY WAVE 


Figure 2.11 Three views of the ordinary and extraordinary reflection points and the phase front of 
a gravity wave. 
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3. DATA COLLECTION AND PROCESSING 

The digital ionosonde, used to provide virtual height data for this 
work is formed by connecting a J5,W ionosonde to a PDP-15 computer through 
a Hewlett-Packard Model 5610A analog-to-digital converter (ADC) . A block 
diagram of the digital ionosonde is shown in Figure 3.1. 

This chapter describes the operation of the system. The next section 
describes the basic hardware. Section 3.2 describes the algorithm used to 
provide accurate virtual height data and the accuracy of the algorithm is 
discussed in Section 3.3. Sections 3.4 and 3.5 discuss the PDP-15 programs. 
Section 3.6 describes a digital frequency counter for the ionosonde and the 
ways it is used. 

3.1 Hardware Interconnections and Timing 

The ionosonde contains a 20 kw pulse transmitter, a receiver and oscil- 
loscopes and cameras to record virtual height data. A pulse generator 
synchronizes the transmitter pulse to the 60 Hz power line and is used to 
synchronize the ADC to the ionosonde. 

The ionosonde transmits a 70 psec pulse every 16.6 msec. The fre- 
quency of the pulse is adjustable from 250 kHz to 20.25 MHz. The bandpass 
of the ionosonde reciever is set by a switch at 10, 20, 35, or 70 kHz. 

The 20 kHz setting, which gives a received pulse width of approximately 
30 psec at 3 dB points, is normally used. The demodulated or video output 
is connected to the ADC through a line driver and a coaxial cable. The 
ionosonde sync pulse triggers the encode gate generator which generates 
the encode gate signal. 300 psec after a pulse is transmitted, the encode 
gate line goes to a logical 1(5 V) for 6 msec as shown in Figure 3.2. 

While the encode gate signal is at "l” the ADC, if enabled by the computer. 



J-5W IONOSONDE 


PDP- 15/30 



I I 


Figure 3.1 Block diagram of the digital ionosonde. 














Figure 3.2 Timing diagram for the digital ionosonde. 
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samples the video output of the ionosonde at 10 ysec intervals, converts the 
analog voltage to a binary value and stores the values in consecutive loca- 
tions in the PDP-15 core memory. The series of samples of the receiver 
output are called a frame of data, 

A program on the PDP-15 identifies the echoes and prints their height 
and amplitude. The basic techniques for identifying echoes and measuring 
their virtual height are discussed in the next section. The data collection 
programs are discussed in Section 3.4. 

3,2 Data Processing Techniques 

A simple method of identifying echoes and determining their height is 
to set a threshold and to call a series of samples greater than the threshold 
an echo. The location of the largest sample in the echo determines the time 
delay of the echo and its virtual height. If the Jth sample in a frame is 
the peak of an echo, the time delay is T = 300 + 10 • (J-l) ysec. 

There are several problems with this method. Noise sometimes exceeds 
the echo amplitude and, depending on the threshold level, either noise is 
detected as an echo or the echo is not detected. To reduce the effects of 
noise, a number of frames are averaged together. 

The received amplitude at a fixed time delay after the transmission of 
a pulse contains both a noise and echo signal. Over a short period of time, 
the signal amplitude is constant but the noise varies in a random manner. 

If the signal exceeds the noise, the average amplitude over a number of 
frames approaches the signal amplitude. The signal- to- noise ratio is in- 
creased by a factor of the square root of the number of frames averaged. 

If too many frames are averaged, however, rapid changes in virtual height 
are not detected. As a compromise, 600 frames are averaged, which reduces 
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the noise by a factor of /600 or 27.6 dB and provides a determination of 
the virtual height every 10 sec. It is preferable to measure the average 
power or average amplitude squared rather than the average amplitude but 
the PDP-15 is not fast enough to square each sample and still have time 
to detect echoes and print out data. 

The resolution to which the time delay is measured is + 5 usee which 
gives a height resolution of only *' 750* m. To increase height resolu- 
tion, the shape of the received echo is interpolated from its samples and 
the center is determined more accurately. Figure 3.3 shows a typical re- 
ceived echo shape and the 10 usee sampling points. A parabola is used to 
approximate the shape of the peak of the received echo which is approxi- 
mately Gaussian. Using the coordinate system in Figure 3.3, the formula 
for a parabola is Y - Ax + Bx + C . The parabola is fitted to the largest 
sample and the sample on each side of the largest. 


at X * 0, Y = y 2 = C 

at 1= 1 ,Y = I 3 = A+ B + C 

at*=-l,Y = Y 1 = 4- B + C 




The maximum amplitude of the parabola is offset from x - 0 by a dis- 
tance P and at this point the derivative of the equation for the parabola 
is 0. 


2 A P + B = 0 



i th SAMPLE 




es of the received echo 


31 


P “ 2A + 2(Y 1 -2I 2 +Y 3 ) 

The echo delay can now be determined more accurately as 

T = 300 + lOfl-l+P) usee 

The accuracy of this formula depends on how close a parabola approxi- 
mates the received echo shape and how much noise distorts the shape of the 
echo. This is discussed in the next section. 

3.3 Error Analysis 

The following analysis of the accuracy of the height algorithm assumes 
that if there are more than one echo, they do not overlap. 

If the center of the echo is offset from the largest sample by a dis- 
tance of Gy the error in determining the virtual height of the echo is P-G. 
Figure 3.4 shows the error as a function of the echo offset G for various 
echo widths. For the 30 usee wide echo, the maximum error is .125 ysec cor- 
responding to a virtual height error of 18.75 m. The error decreases with 
increasing pulse width but, because the wider the echo, the closer the 
three samples are to being the same amplitude, noise causes larger errors 
for wider echoes than for narrow echoes . 

To simulate the error under different noise conditions, a Monte Carlo 
approach is used. Random noise with a Gaussian distribution is added to a 
Gaussian shaped echo. A parabola is fitted to the largest sample and the 
sample on either side and the offset P is calculated. The process is re- 
peated many times to get an RMS error. Figure 3.5 shows how the error 
varies with offset for a 30 ysec wide echo under conditions of no noise, 

50 dB and 30 dB signal-to-noise ratio. The error for less than 50 dB 
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Figure 3,4 Height error as a function of G for various echo widths and no noise. 
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Figure 3.5 Height error as a function of G for a 30 ysec wide echo and various 
noise conditions. 
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signal -to-noise ratio is approximately independent of the offset G . Figure 
3.6 shows the error as a function of the received pulse width for an offset 
of 2.5 ysec and various signal conditions. For a fixed signal-to-noise 
ratio, the error has a minimum value for some range of pulse width. 

Noise conditions for the ionosonde vary with frequency, time of day 
and weather. Typical signal-to-noise ratios are from 0 to 30 dB. Averaging 
over 600 samples adds about 30 dB to the signal-to-noise ratio. Figure 3.6 
indicates that the 30 ysec pulse width is near the minimum error for signal- 
to-noise ratios in the range of 30 to 60 dB. The error varies from + 300 m 
in noisy conditions to + 50 m for good signal conditions. 

The phase of the returning echo relative to the transmitted signal and 
the time delay for the echo depend on the height of the reflector. When two 
echoes are received from different parts of the reflector but with nearly 
equal virtual heights, the amplitudes of the two echoes add vectorial ly and 
both echo shapes are distorted slightly. The apparent height and amplitude 
of each echo are different from the actual height and amplitude if the 
echoes did not interfere. The amount of change depends on the phase dif- 
ference between the two echoes, their separation and their relative ampli- 
tudes. The maximum changes occur when the two echoes are in phase and 180° 
out of phase. 

For the lower of the two echoes, the apparent height and amplitude 
increases when they are in phase and decreases when they are 180° out of 
phase. For the higher echo, the apparent height decreases and the ampli- 
tude increases when they are in phase. The height increases and the ampli- 
tude decreases when they are 180° out of phase. The changes in echo shape 
are shown in Figure 3.7. The maximum height and amplitude changes as a 
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Figure 3.6 Height error as a function of echo width with G=2.5 ysec and various noise conditions. 
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Figure 3,7 Distortion of echo shape due to two overlapping echoes. 
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function of the separation between the echoes for several values of rela- 
tive amplitude are shown in Figure 3.8. 

3.4 Collection Programs 

A set of programs RPROC, HDATA, ECHO and PCK are used to collect and 
print out virtual height data on the digital ionosonde. An overall flow 
chart of the collection process is given in Figure 3.9. 

RPROC is the main program which calls other subroutines. HDATA is a 
collection of subroutines which are written in assembly language because 
their functions cannot be done in FORTRAN. HDATA contains the following 
subroutines. 

Subroutine GO reads console switch 0 which determines whether collec- 
tion continues or restarts. 

Subroutine ADGET sets up the ADC to put 530 samples in array BUF1. 

When the frame of data is in BUF1, an interrupt generated by the ADC starts 
the interrupt service routine ADINT at the highest priority level. Figure 
3.10 contains a flow chart of the interrupt service routine. ADINT checks 
for a timing error in the ADC and starts the real-time subroutine RTS at a 
lower priority level. Additional operations done by RTS are not done in 
ADINT because the disks need to interrupt on the highest priority level 
but cannot while ADINT is running. 

If a timing error occurs in the ADC or if the first sample in BUF1 is 
above a fixed level indicating a noisy frame, the data on that frame is 
ignored, otherwise the elements in BUF1 are added to the elements in array 
BUF2. After 600 frames, minus those ignored, are added together in this 
manner, the data in BUF2 is transferred to array CDATA and BUF2 is zeroed 
out to prepare for averaging the next 600 frames. The code used to add the 
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Figure 3.8 Height and amplitude change as a function of echo separation for several 
values of relative amplitude. 
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Figure 3.9 Overall flowchart of data collection and processing 
programs . 
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Figure 3.10 Flowchart of the interrupt service routine. 
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600th frame to BUF2 and transfer BUF2 to CDATA is executed in less than 
11.4 msec in order to start the ADC in time to catch the next frame of data 
as shown in Figure 3.2. The code to add BUF1 to BUF2 is in straight line 
form to minimize its execution time. Since this code is executed once for 
each frame, minimizing its execution time provides more time for processing. 

Subroutine RDATA is used to transfer control to the background until 
600 frames have been added together, then data processing by RPROC continues. 
A flowchart of subroutine RDATA is given in Figure 3.11. Occasionally the 
ADC does not start when requested to by ADGET. RDATA checks for this condi- 
tion and restarts the collection process if it has stopped. If the data 
processing by RPROC takes longer than 10 sec, the data collection and pro- 
cessing get out of synchronization. RDATA checks for this condition also 
and prints a warning if it occurs. 

Subroutine RDAT sets up the inline code used to add BUF1 to BUF2, zeros 
out BUF2 and starts the ADC at the beginning of a collection run. 

Subroutine BGDAT sets up a series of words which contain pointers to 
data used by the background display programs discussed in Section 3.5. 

The remaining programs, PROC, ECHO and PCK, are written in FORTRAN, 

Subroutine ECHO searches a copy of part of CDATA for echoes. A flow- 
chart of ECHO is given in Figure 3.12. ECHO searches a region from 75 km 
below to 75 km above a specified height. The region searched is wide 
enough to include the ordinary and extraordinary echoes under most condi- 
tions. An echo is defined as 4 or more continuous samples above a threshold 
and noise is defined as all samples below the threshold. A threshold is 
set initially and a new threshold is computed for each 10 sec of data. 

Because of the offset in the amplifiers between the ionosonde and the ADC, 
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Figure 3,11 Flowchart of subroutine RDATA. 
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Figure 3.12 Flowchart of subroutine ECHO. 
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the true zero level is not known. The sample with the smallest value is 
assumed to be the zero level. The average of the noise samples is calcu- 
lated and a new threshold is set at 5 times the average noise level. The 
threshold adjusts to changing noise and signal conditions. If the noise 
level suddenly increases above the threshold between one 10 sec period and 
the next, when the average noise is calculated, no samples will be below 
the threshold and an improper noise level is obtained. If too few samples 
are below the threshold, the threshold is set at a high level and another 
attempt is made to determine the noise and compute a new threshold. 

Once an echo has been found, a parabola is fitted and the amplitude 
and height are recorded. Then a modified parabola shape shown in Figure 
3.13 is subtracted from the actual echo shape and the search starts again 
below the location where the last echo started. This allows echoes which 
are close together to be resolved. 

When RPROC is started, an initial height XHT is typed in. ECHO 
searches for echoes near this height and again at twice this height for 
second reflections. RPROC prints out data on 3 echoes at the lower height 
and 2 echoes at the higher height. The data are also put in a compact 
form and stored on the disk. After the echoes near XHT are located, sub- 
routine PCK is called. 

Subroutine PCK looks for an echo within 20 km of XHT. If there is at 
least one, the height of the closest echo to XHT is used as the new value of 
XHT. This is done so that ECHO is always searching a region centered on 
the location of the last echo. 

3.5 Display Programs 

In February 1973, a Hazeltine 2000 cathode-ray terminal replaced the 
model 33 teletype as the background control device. The Hazeltine terminal 
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Figure 3.13 Subtract assumed echo shape from received echo to resolve second echo. 
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operates on a serial current loop interface and is compatible with the 
interface for a teletype. The Hazeltine terminal operates at a maximum 
speed of 9600 baud or 873 characters per second as compared to the tele- 
type which operates at 110 baud or 10 characters per second. The sub- 
stitution required only minor modifications to the PDP-1S hardware to use 
the full speed of the Hazeltine terminal. 

The background clock board M452 is shown in Figure 3.14. The oscil- 
lator generates an 800 Hz signal which is divided to provide 110 Hz. 

Capacitor C5 is decreased to .01 pf and resistor R3 is adjusted to give 
76800 Hz for the required baud rate of 9600. The TTY connector W078A does 
not operate at this rate because of a 1 pf filter capacitor. After changing 
this capacitor to .022 yf the Hazeltine terminal worked properly. 

The Hazeltine terminal writes 5 full 74 character lines per second. 

The entire screen of 26 lines can be filled in approximately 5 sec. The 
Hazeltine terminal is used to provide real-time information on the opera- 
tion of the data collection process and on the data being collected. 

Figure 3.15 contains a flowchart of the background display programs BGRAPH, 
BDATA, PARA and GRAPH. Depending on the setting of console switch 17, either 
a graph of 11 minutes of virtual height data or a graph of the echo shape is 
produced. An example of the virtual height graph is shown in Figure 3.16. 

A starting time in minutes is typed in and then the height of the echo 
picked by PCK is displayed. This graph provides a quick method of deter - 
mining what the reflecting layer is doing and whether collection should 
continue or not . 

An example of the echo-shape graph is given in Figure 3.17. The 
height and time for the one-hop echo in the center of the screen is given 
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Figure 3.14 Schematic diagram of background clock board M452 
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Figure 3.16 Example of height versus time graph. 
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at the top on the graph. The 0's represent the one-hop echo, the +'s repre- 
sent a parabola fit to the one-hop echo and the #’s represent the two-hop 
echo. The left-hand side of the screen has the height and amplitude of the 
averaged data near the one-hop echo. If any of the symbols overlap, an * 
is printed. The one-hop threshold is indicated by a T on the third line of 
the graph. This graph is useful in determining whether the algorithms in 
ECHO are functioning properly. It is also useful to check that ECHO is de- 
tecting all the echoes and not detecting noise as echoes. In conditions of 
high noise and low- signal strength, echoes can be seen on the graph that 
cannot be seen on an oscilloscope monitoring the video signal from the 
ionosonde. 

3.6 Ionosonde Frequency Counter 

A frequency counter was designed and built to measure the ionosonde 
transmitter frequency. The ionosonde is licensed to operate only on the 
ten frequency bands listed in Table 3.1. The frequency counter keeps the 
operation of the ionosonde legal by turning the transmitter off while the 
ionosonde is sweeping through frequencies not listed in Table 3.1. A digi- 
tal display of the ionosonde frequency is provided in the computer room 
and, with remote ionosonde controls, the ionosonde can be operated from 
the computer room. With the digital input port to the PDP-15, the com- 
puter reads the ionosonde frequency and by using the digital output port, 
the computer can control the ionosonde frequency. 

The ionosonde frequency range, 250 KHz to 20.25 MHz, is divided into 
four bands, 250 kHz to 750 kHz, 750 kHz to 2.25 MHz, 2.25 MHz to 6.75 MHz 
and 6.75 MHz to 20.25 MHz, called Band I, II, III and IV respectively. 

The transmitter frequency is formed in a mixer which subtracts a 4 MHz 
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signal on Bands I, II and IV or 8 ttlz on Band III from a variable fre- 
quency oscillator (VFO) signal. The mixer operates in a pulsed mode so 
the transmitter frequency is not continuously available and cannot be 
counted directly. The counter measures the VFO frequency minus 4 or 8 MHz. 
This is done by presetting the counter to - 4 MHZ (- 8 MHz on Band III) 
before counting starts . 

The circuit diagram for the frequency counter is given in Figures 3.18, 
3.19 and 3.20. A 1 MHz crystal oscillator is divided by 400 to provide the 
time base. The WO signal is converted to TTL levels and counted for a 
period of 200 ysec which provides a frequency resolution of 5 kHz. The 
counter is 12 bits long, giving a maximum frequency count of 20.48 MHz. 

The timing diagram in Figure 3.21 shows the relation between the preset 
period and the counting period. During the preset period the counter is 
preset to 6340 0 (4700„ on Band III) and the SRFF formed by two cross -coupled 
NANDs is set. 

If the VFO input signal is not present, during the clear period the 
counter is set to 0. If this were not done, when the ionosonde is off, the 
counter would indicate 6340 g . When the computer reads a frequency of 6340 g , 
it cannot tell whether the ionosonde is operating at 16.48 MHz or turned 
off. A reading of 0 is unambiguous. If a VFO signal is present, the 
counter will count during the count period, the SRFF is reset and the 
counter is not cleared during the clear period. At the end of the count 
period, the counter is allowed to settle for 50 ysec, then the Data-Ready 
signal is true for 50 ysec. During this time, the count is transferred to 
12 latches. The data in the latches are read by the computer through the 
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Schematic of ionosonde frequency counter time base, 


counter and frequency deleter board. 
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Figure 3.19 Schematic of ionosonde frequency counter data latch and binary to BCD converter board. 
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Figure 3.20 Schematic of ionosonde frequency counter power supply and board interconnections 
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The data from the latches are converted to 4 decades of BCD using 
7 code converter read-only-memories (ROMs) . The BCD data drive a frequency 
display in the computer room, indicating the ionosonde frequency in 10' s 
of kHz . 

The 10 most significant bits of the counter are used to determine if 
the ionosonde frequency is in an allowed frequency band or not. The fre- 
quency range from 0 to 20.48 MHz is covered in 1024 segments, each 20 kHz 
wide. A 1024 bit ROM has one bit for each segment. A "l" is programmed 
in the ROM for segments that are in prohibited frequency bands and a "0" 
is in the other bits. 

The 1024 bit ROM used in the counter is organized as 256 four-bit words. 
A 4-bit data selector effectively converts the ROM into a 1024 by 1 organiza- 
tion. If the output of the data selector is 1 during the Data-Ready period, 
a monostable multivibrator is triggered and puts out a 33 msec pulse used 
to turn off the ionosonde transmitter. If the ionosonde is in a restricted 
frequency, the monostable multivibrator is retriggered every 400 ysec and 
stays set, keeping the ionosonde transmitter off. If the ionosonde fre- 
quency is on the edge of an allowed band, because of the uncertainty of + 1 
in the count, the monostable multivibrator is not retriggered every 400 ysec 
but once every several 400 ysec periods. As long as the monostable multi- 
vibrator is triggered at least once every 33 msec, it stays on and keeps the 
transmitter off. 

The monostable multivibrator Q and Q outputs drive 2 high current inter- 
face drivers to light LED's, indicating whether the transmitter is on or off 
and also drive a line which turns the ionosonde transmitter on and off. A 
switch is provided to override the counter and keep the ionosonde constantly 
on or off. 
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4. EXPERIMENTAL RESULTS 

This chapter presents an analysis of some of the data collected on 
the digital ionosonde by methods described in Chapter 2. 

The digital ionosonde was operated for about two hours starting 
between 1300 and 1500 CST several days a week for five months. A fre- 
quency of 4.1 MHz was normally used because there is little interference 

/ 

at that frequency and good reflections from the F region are obtained. 

The extraordinary reflection cannot be seen at this frequency and 
several runs were made at a higher frequency where both reflections were 
present. 

On each run, the virtual height varies in a quasiperiodic manner 
with a period of 10 min or more and a height change of more than 10 km 
and in some cases as much as 50 km. Due to the increasing height of the 
F region in the afternoon, the fluctuations are superposed on a linear 
increase of virtual height. The amplitude generally fades with a period 
of a few seconds because of very small irregularities in the ionosphere. 
This fading is nearly all removed by the present system which averages 
over a 10 sec period. Their residual effect is to produce small random 
height and amplitude variations. Some of the random variation in height 
is due to errors in the height algorithm described in sections 3.2 and 
3.3. Large periodic increases in amplitude occur in some runs with periods 
of 15 to 30 min due to focusing. Focusing is observed on 75% of the data 
runs and several of these are described in section 4.3. 

The multiple reflection points described in section 2.4 are shown in 
33% of the runs. Section 4.4 shows examples of these. 

Data taken at frequencies higher than 4.1 MHz, showing both ordinary 
and extraordinary echoes are described in section 4.5. 
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Section 4.6 describes some experimental data with features that can- 
not be explained with the present reflector model described in Chapter 2. 

In the next section, a sinusoidal reflector is used to illustrate 
the effects of the length of the reflector on focusing and multipath. 

Section 4.2 describes a data run in which two echoes are close 
enough to interfere with each other as described in section 3.3. 

4.1 Effects of Reflector Horizontal Wavelength on Focusing and 
Multipath 

Focusing is caused by the reflector acting as a curved mirror as 
described in section 2.3. Multipath is when more than one echo is 
observed from one reflector simultaneously as described in Section 2.4. 

To illustrate the effects of the reflector wavelength on the occurrence of 
focusing and multipath, the virtual height and amplitude for sinusoidal 
reflectors of various wavelengths are calculated The reflectors have 
the form g - 250 - 10cos(2ttx/P) km. The average height and amplitude of 
the reflectors are typical values for a frequency of 4.1 MHz. P is the 
wavelength of the reflector. The amplitude and virtual height of an echo 
from the reflector are calculated for one cycle of the reflector, 0 < x < P. 

The horizontal wavelength of the reflector in Figure 4.1 is 500 km. 
Only a small amount of focusing occurs in the two-hop amplitude and almost 
none in the one-hop amplitude. The virtual height of the one-hop comes to 
a narrower peak than the reflector and the two-hop height is more sharply 
peaked than the one-hop height. For longer wavelengths, there is less 
focusing and the virtual height is closer to the reflector height. Near 
the center of the reflector, the two-hop echo experiences more focusing 
than the one-hop echo for two reasons. First, the two-hop echo is re- 
flected and focused twice, effectively doubling the depth of the reflector 
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Figure 4*1 One- and two-hop virtual height and echo amplitude for a sinusoidal reflector 
with P - 500 km. 
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for the two-hop echo. Secondly, the two-hop echo has to travel twice as 
far making the apparent reflector for the two-hop echo twice as high. For 
long wavelengths, the transmitter and receiver are closer to the center 
of curvature for the apparent two-hop reflector than the one-hop reflector. 

As the reflector wavelength decreases to 450 km, the two-hop height 
comes to a sharp peak and the amplitude shows a large enhancement as 
shown in Figure 4.2. The one-hop height and amplitude are not much dif- 
ferent from the height and amplitude in Figure 4.1. 

For a wavelength of 400 km, as shown in Figure 4.5, three reflec- 
tion points appear in the two-hop virtual height. The one-hop height 
shows increased distortion and the amplitude shows some focusing. 

For a wavelength of 300 km, as shown in Figure 4.4, three reflec- 
tion points appear on the one-hop virtual height. 

Except for the cases of multipath, the maximum and minimum one-hop 
virtual height correspond to the maximum and minimum height of the 
reflector. 

4.2 Interference Between Two Echoes with Nearly the Same Virtual Height 

The interference between two echoes is described theoretically in 
Section 3.3. This section describes an example of this phenomenon. 

Figure 4.5 is a graph of 100 min of data taken from 1931 to 2111 CST 
on May 4, 1973 at 3.1 MHz with almost four cycles of a disturbance with a 
period of 26 min. The two-hop echo shows four segments of focusing but 
at other times the two-hop amplitude is too small to be detected. For the 
first 50 min, the one-hop height has a random fluctuation of about 160 m 
superimposed on the slower periodic variation. This error is within the 
accuracy range of the echo detection algorithm as calculated in Section 3.3. 
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Figure 4.2 One- and two-hop virtual height and amplitude for a sinusoidal reflector with P = 450 km. 
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Figure 4.3 One- and two-hop virtual height and amplitude for a sinusoidal reflector with P = 400 km. 
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Figure 4,4 One- and two-hop virtual height and amplitude for a sinusoidal reflection with P = 300 km. 
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An interesting feature of this record is the short -period oscillation 
in the one-hop height and amplitude from 2022 to 2039. This portion of the 
graph is expanded in Figure 4.6. The same feature is noted again from 
2044 to 2111, but not as clearly. While the data were being taken, two 
echoes very close together were noted visually, but the echo detection 
algorithm described in section 3.2 which can detect two echoes close 
together was not yet in use. The detection algorithm used finds only one 
echo which is the echo within 10 km of the echo found in the previous 
10 sec interval if such an echo exists. If more than one echo is within 
10 km of its previous echo, the closest echo to the echo found in the pre- 
vious 10 sec is picked and printed out . Only the lower echo was detected 
except at 2023 and 20 sec when the upper echo was detected. The difference 
in the virtual height of the two echoes decreases from 2022 to 2030 and 
then increases . 

The phase difference between the echoes changes from 0° to 180° and 
back to 0° every time the difference between their virtual heights changes 
by one-half wavelength. When the echoes are in phase, the apparent height 
and amplitude of the lower echo increases and when they are out of phase, 
the amplitude and apparent height decreases. At 2030, the two echoes are 
closest together. One minute later, the apparent height and amplitude 
have gone through one oscillation and the echoes have moved one-half wave- 
length apart or 48.4 m. The frequency of the oscillations increases indi- 
cating that the rate of separation is increasing. At 2035, the oscillations 
appear to die out. Because the height and amplitude are averaged over 
10 sec, what is actually happening is that the height and amplitude are 
going through one complete oscillation in 10 sec. The rate of separation 
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Figure 4.6 Expanded view of part of Figure 4.5 showing interference between two echoes 
and the calculated position of the higher echo. 
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continues to increase until the echoes are too far apart to interfere. The 
same thing appears to happen while the echoes are approaching but it is 
not as clearly represented in the graph. 

From the amount of change in amplitude and height and using Figure 3.8, 
the approximate virtual height of the second echo near 2030, where the 
oscillations are slow enough that the 10 sec averaging does not distort the 
oscillations, is determined and put on Figure 4.6. The number of phase 
variations determines the height of the upper echo at other times. The 
calculated position comes close to the one point where the upper echo was 
detected. The height of this point is also distorted by as much as the 
lower echo height. The detected point is within this error of the calcu- 
lated echo position. 

The difference between the radial velocities of the two echoes is 
very small, on the order of a few m sec - " 1 , indicating that their motion is 
coupled in some way. A possible explanation for the two echoes is that 
the shape of the reflector is not constant in the Y direction as assumed 
in section 2.3. Another possibility is that a small additional irregu- 
larity exists at the same height as the reflector but displaced several 
km laterally from the ionosonde antennas in the direction of propagation 
of the disturbance. The irregularity is displaced vertically by the dis- 
turbance and follows the height of the reflector. As the reflector slope 
increases, the reflection point moves in the direction of propagation of 
the disturbance and comes closer to the irregularity until, at the inflec- 
tion point of the reflector, halfway between the minimum and maximum of 
the reflector, the reflection point is at its closest approach of the 
irregularity. This is the only occurrence of this phenomenon in five 
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months of data collection but the same effect occurs when multipath gives 
two echoes close together. 

4»3 Data Shouting Focusing 

Focusing, as described in Section 2.3, is observed on three-quarters 
of the data runs taken. Several of these are described in this section. 

Figure 4.7 shows 65 min of data taken on May 5, 1973 at 1502 CST. 

Two cycles of a disturbance are shown with a period of about 22 min. The 
amplitude shows large enhancements just after each peak in height. The 
two-hop echo is not included in Figure 4.7 because it is only present 
twice when focusing increases the two-hop echo amplitude to a point where 
it is detected. Using the calculator program described in Section 2.4, 
the reflector shape that produces the focusing effects in Figure 4.7 is 
determined as shown in Figure 4.8. If the overall dimension of the re- 
flector is reduced to 1700 km, multiple reflection points form. If it is 
increased to 1900 km, the focusing produced by the first cycle disappears. 
The dimension of the reflector must be somewhere between 1700 and 1900 km. 
1800 km provides a good fit to the data. The wavelength of the distur- 
bance is then 900 km and the velocity is 680 m sec ^ . 

Figure 4.9 is a graph of 83 min of data collected on November 2, 1973 
at 1336 CST. Three cycles are shown with an average period of 24 min. 

The extraordinary echo is detected at several points during the run. The 
amplitude for the one-hop echo shows rapid fading but no clear focusing. 

The height has a random variation of 630 m superimposed on the periodic 
variations. Small irregularities in the ionosphere are causing some of 
the random height variation because the height accuracy of the digital 
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Figure 4.8 Reflector shape that produces virtual height and amplitude curves that match 
the data in Figure 4.7. 
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Figure 4.9 Data record taken on Nov. 2 , 1972 showing focusing of the 
two-hop echo. 
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ionosonde is greater than this. Signal conditions were very good for this 
run and the two-hop echo is detected for most of the run. The two-hop 
amplitude shows two periods of focusing. Figure 4.10 shows the reflector 
shape and the one- and two-hop height and amplitude calculated from the 
reflector. A reflector 1000 km long produces the best fit to the data. 

This gives a wavelength of 475 km and a speed to 396 m sec -1 . 

4.4 Data Showing Multipath 

Multipath occurs when more than one one-hop or two-hop echo is 
received simultaneously from the reflector as described in Section 2.4. 

This section presents and discusses several data runs showing this effect. 

Figure 4.11 shows 122 min of data collected on September 7, 1973 at 
1308 CST and at a frequency of 4.1 MHz. There appear to be 13 cycles of 
a disturbance in the first 90 min giving an average period of 7 min. The 
disturbance appears to die out and be replaced by a slower disturbance. 

Two cycles at 1344 and 1354 show multipath. These are expanded in Figure 
4.12. From 1344 to 1345 the two echoes near 240 km have nearly the same 
virtual height and are close enough to interfere with each other as des- 
cribed in Section 3.3. The echo heights in Figure 4.12 are distorted 
because the two echoes are close together. The same thing happens with 
the two echoes near 250 km from 1354 to 1356. The reflector which produces 
these echoes is shown in Figure 4.13. The wavelength of the disturbance 
is about 725 km which gives a velocity of 450 m sec * 

Figure 4*14 show 135 min of data taken on December 8, 1973 starting 
at 1309 CST and at a frequency of 3.8 MHz* About 6 cycles of a disturbance 
with a period of 15 min are shown* The ordinary and extraordinary echoes 
are both recorded and also a number of two-hop echoes* The one-hop echoes 
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Figure 4.10 Reflector shape that produces virtual height and 

amplitude curves that match the data in Figure 4.9 
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Figure 4.13 Reflector shape that produces virtual height and amplitude curves that match 
the data in Figure 4.12. 
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show little focusing and no splitting. The two-hop echoes are split in 
at least three places at 1328, 1357 and 1423. The two-hop amplitude is 
generally too low to give a good picture of what is happening. Figure 4.15 
shows a reflector shape which produces two-hop echoes similar to those in 
Figure 4.14. From this reflector, the approximate wavelength of the dis- 
turbance is determined to be 400 km and the velocity of the disturbance 
is 456 m sec ^ . 

4.5 Data Shouting Ordinary and Extraordinary Reflections 

This section shows how the ordinary and extraordinary echoes can be 
used together to determine the direction of the disturbance. Two data 
records with both ordinary and extraordinary echoes present are analyzed 
in this section. The basis for the direction measurement is explained in 
Section 2.6. 

Figure 4.14, used as an example of multipath, has both ordinary and 
extraordinary echoes. The extraordinary reflector is delayed by approxi- 
mately 55 sec from the ordinary reflector which means the disturbance is 
traveling roughly in a southerly direction. The frequency used for this 
run was 4.8 MHz and the F region critical frequency was 7.5 MHz. Using 
the program to compute lateral deviation described in Appendix IV, the 
lateral deviation is computed as 9 km and the vertical separation of the 
reflection points is 16.5 km. The period of the disturbance is approxi- 
mately 10 min so a value of 45 deg is assumed for the angle of the phase 
front of the disturbance relative to vertical. Using the formula for the 
direction of the disturbance derived in Section 2.6 and the velocity of 
456 m sec * calculated in Section 4.4, the angle of propagation of the 
disturbance relative to north is 19 deg. The angle of propagation is very 
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Figure 4.15 Reflector shape that produces virtual height and 

amplitude curves that match the data in Figure 4.14 






sensitive to changes in the angle of the phase front. Decreasing the 
angle to 43.7 deg changes the direction of propagation to directly south. 

Figure 4.9, used as an example of focusing in Section 4.3, also shows 
parts of the extraordinary reflector. The variations in height of the 
extraordinary reflector are delayed by about 166 sec on the average from 
the ordinary height variations. This means the disturbance is traveling 
roughly from north to south. The frequency used for this run is 4.2 MHz 
and the critical frequency is 9.5 MHz, Using the program to compute 
lateral deviation, the lateral deviation is 4.12 km and the vertical sepa- 
ration of the reflection points is 8 km. The period of the disturbance 
is about 25 min. Using a phase front angle of 45 deg, and assuming the 
disturbance is traveling directly to the south, the projection of the 
velocity on the north-south axis is 73 m sec ^ which is much slower than 
the velocity of the disturbance calculated as 396 m sec”^ in section 4.3. 
The projection of the velocity on the north-south axis should be larger 
than the magnitude of the velocity or equal to it if the disturbance is 
traveling due south. A value of 82 deg for the phase angle is required 
to make the projection of the velocity equal the magnitude of the velocity. 
This is very close to the angle of 83 deg measured by Georges [1967] for 
a disturbance with a period of 40 min. 

4.6 Data Showing Unusual Features 

Several data runs contain features which cannot be explained with 
the present reflector model described in Chapter 2. 

Figure 4.16 is a graph of 125 min of data taken on November 17, 1973 
at 1340 CST at 4.8 MHz. The amplitude of the one-hop echo is very high 
and during the first 60 min of the run, the ADC sometimes is saturated. 
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An ordinary and extraordinary reflection are present and the ordinary two- 
hop echo is also observed. About 11 cycles of a disturbance with a period 
of 10 min is observed. The amplitude of the disturbance is below normal 
and because of the saturation of the ADC, the shape of the disturbance is 
not clear. The two-hop echo shows segments of focusing which correspond 
to peaks in the ordinary one-hop height but for the most part, the ampli- 
tude is too small to get a good picture of the two-hop height variation. 

An interesting feature of the two-hop height is the echo which in- 
creases linearly in height. The amplitude of this echo increases slowly 
until about 1440 and then decreases. The amplitude is very steady, showing 
no signs of fading. The echo is first picked up at 1359 and 20 sec at a 
height of 397.5 km and increases in height beyond the detection range of 
ECHO at 1502 and a height of 535 km. The low height of the echo from 1400 
to 1430 implies that it is not a two-hop echo since there is no one-hop 
echo at one-half this height. The amplitude is very constant and higher 
than the two-hop echo amplitude when it is not focused. This signal could 
be a type of noise reported by Gerson [1962] which he calls 8 type sweepers. 
The signal could be caused by another pulse transmitter operating at nearly 
the same frequency and synchronized to a power line with a slightly dif- 
ferent frequency to account for the changing time delay relative to the 
ionosonde transmitter. 

Figure 4.17 shows 93 min of data taken on December 1, 1973 starting 
at 1310 CST at a frequency of 4.1 MHz. A few cycles of a disturbance with 
a period of about 20 min are shown. A few extraordinary echoes are also 
detected. The two-hop echo is present for most of the run. An unusual 
feature of this record is the second two-hop trace starting at 1417 and 
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Figure 4.17 Data record taken on Dec. 1, 1973 showing an unusual 
two -hop echo. 
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increasing in height until the end of the run. This trace is similar to 
the one in Figure 4.16 but it is not as straight and does not start as far 
below the normal two-hop height. The amplitude of the echo is approximately 
constant for the 25 min that is recorded. It also starts too low to be a 
two-hop echo. This echo appears to be associated with one cycle of the 
disturbance and it could be that this is a multipath phenomenon which can- 
not be explained with the present reflector model. 

Figure 4.18 shows 122 min of data taken on October 10, 1973 starting 
at 1250 CST at a frequency of 4.7 MHz. This graph clearly shows several 
examples of multipath. There is very little fading in the amplitude and 
the random height error is less than 200 m. Focusing can clearly be 
seen on each cycle of the disturbance. The two-hop echo also shows 
focusing but for the most part the two-hop echo is too weak to be de- 
tected. A few extraordinary echoes are detected. At this frequency and 
on the day the data were taken, the extraordinary echo is below the 
ordinary echo. The exact shape of the multipath echoes cannot be recon- 
structed with the present curved reflector model. All multipath echoes 
using the current reflector model have some form of X or crossover which 
is not present in Figure 4.18. 
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5. RESULTS AND SUGGESTIONS FOR FUTURE WORK 


In the previous chapters, a system for collecting accurate fixed- 
frequency virtual-height data and methods for analyzing the data to 
determine speed and direction of traveling ionospheric disturbances are 
presented. Virtual height is measured with a time resolution of 10 sec 
and with a height resolution of less than 300 m and in some cases as 
little as 50 m. Five months of data from the digital ionosonde shows 
quasi-periodic height variations on all data records with periods of 5 
to 30 min and amplitudes of 10 to 50 km. Focusing is apparent on 75% 
of the records and multipath is observed on 33% of the records. 

The horizontal wavelength of the disturbance is determined by 
matching virtual height and amplitude calculated from an assumed re- 
flector shape with the data from the digital ionosonde system. Several 
data records analyzed in this manner indicate wavelengths from 400 to 
900 km and speeds of 400 to 680 m sec -1 . The direction of the distur- 
bance is determined by measuring the time delay between the appearance 
of prominent features of the disturbance at the ordinary and extra- 
ordinary reflection points. Two records analyzed in this manner indicate 
that the disturbance is traveling approximately due south. 

To graph the data from the digital ionosonde on the plotter, large 
amounts of paper tape are punched out on the PDP-15 and generally used 
only once to make a graph of the data. The partial reflection experiment 
and presumably future experiments will collect data on the PDP-15 and 
graph it on the plotter. The paper punch on the PDP-15 occasionally makes 
errors in punching the tape and the paper tape reader connected to the 
calculator also makes occasional errors in reading the tape. The tape 
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reader also has no convenient way of handling large pieces of paper tape. 
For these reasons, another method of transferring data between the PDP-15 
and the calculator is desirable. A magnetic cassette tape unit could be 
attached to the PDP-15 through the digital input-output port. This would 
eliminate the need to use paper tape for transferring large amounts of 
data. It would be necessary to find out the format used by the calculator 
in writing on the cassette tapes and program the PDP-15 to write tapes 
in the same format. 

The frequency counter described in Section 3.6 makes computer control 
of the ionosonde frequency possible. It should be possible to switch the 
ionosonde between several frequencies to get a vertical profile of the dis- 
turbance. The frequency counter also simplifies writing a program to take 
digital ionograms. 

To investigate the unexplained data and increase the accuracy of the 
wavelength measurements, ray tracing could be used with a model of the 
ionosphere as disturbed by gravity waves. Clark [1970] develops a com- 
puter model of the ionosphere disturbed by gravity waves and a set of 
differential equations for ray tracing by Haselgrove [1963] is implemented 
on the HP-9830 calculator. A data collection program that takes continuous 
digital ionograms or one that takes an ionogram every few minutes and 
fixed frequency virtual height data at other times could be used to 
determine the electron density distribution. 

A first step in improving the data from the digital ionosonde is to 
add phase measuring capability by adding two synchronous detectors to the 
ionosonde. The computer, using the digital output port sets the height 
at which two sample-and-hold circuits sample the output of the synchronous 



detectors. After the ADC finishes sampling the magnitude of the ampli- 
tude as it does in the present system, the two sample-and-hold voltages 
are multiplexed into the ADC input and sampled. 

A further improvement is to measure the angle of arrival of the echo 
Two more antennas, two more receivers, four more synchronous detectors 
and four saraple-and-hold units would be required. Since the computer can 
rapidly change the height at which the sample and hold circuits sample, 
if there are more than one echo, a different echo can be looked at with 
each transmitter pulse. The direction of propagation of the disturbance 
is obtained directly and unambiguously from the angle-of-arrival 


measurement . 



91 


REFERENCES 

Clark, R. M. , K. C. Yeh and C. H. Liu (1970), Gravity waves in the iono- 
sphere, Tech. Rep . 39, Ionospheric Radio Laboratory, Univ. 111., 

Urb ana-Champ ai gn . 

Budden, K. G. (1960), Radio waves in the ionosphere, Cambridge University 
Press, London. 

Georges, T„ M. (1967), Ionospheric effects of atmospheric waves, ESSA 
Tech. Rep . IER 57 -ITS A 54. 

Gerson, N. C. and W. H. Gossard (1962), Diurnal variation of sweeper 
activity, J. Geophys. Res. 77, 1007-1011. 

Harper, R. M. (1972), Observation of a large nighttime gravity wave at 
Arecibo, J . Geophys. Res . 77, 1311-1315. 

Haselgrove, J. (1963), The Hamiltonian ray path equations, J. Atmos. Terr. 
Phys., 25, 398-400. 

Hines, C. 0. (1954), An interpretation of certain ionospheric motions in 
terms of atmospheric waves, J. Geophys. Res. 64, 2210-2211. 

Munro, G. H. (1950), Traveling disturbances in the ionosphere, Proc. Roy. 
Soc. 202, 208-223. 

Whitehead, J. D. (1960), Focusing of short radio waves reflected from the 
ionosphere, J. Atmos. Terr. Rhys., 9, 269-275. 

Yeh, K. C. and C. H. Liu (1972), Theory of Ionospheric Waves, International 
Geophysics Series 17, Academic Press, New York. 



92 


APPENDIX I. OPERATING INSTRUCTIONS FOR THE DIGITAL IONOSONDE AND A 
LISTING OF THE DATA COLLECTION AND DISPLAY PROGRAMS 

The following procedure is used to operate the fixed frequency data 
collection system. A listing of the programs is included after the 
operating instructions. 

1. Turn on Hazeltine 2000 video terminal and the PDP-15 computer. 

Press the ONLINE button on the Terminet 1200 TTY. The address 
switches on the PDP-15 should be set to OOlllg. Press STOP and 
RESET at the same time. Then press START. The following should 
appear on the foreground TTY. 

FKM15 V3A 

$} 

If not, see the procedure for cold starting the PDP-15 taped next 
to the computer console. 

2. Turn on the encode gate generator and monitor oscilloscope. 

Disconnect the coaxial cable from the PR receiver to the ADC input 
and connect the cable from the ionosonde video attenuator to the 
ADC input. The cable from the PR encode gate generator is simi- 
larly disconnected from the ADC encode input and replaced with the 
cable marked IONOSONDE ENCODE PULSE. 

3. Mount the tape marked Ionosonde Programs on tape unit 2. At the 
foreground TTY type the underlined portion of the following. 

$} N 2 cr cr means press the CARRIAGE RETURN 

$} A DT2 -4, -5/DK2 2 cr key 

$}MPOFF cr 


$}SHARE cr 
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$ } GLOAD cr 
FGLOAD V2A 

$ } PROC,HDATA esc esc means P ress the ESCAPE key 

At the background console type the following. 

BKM15 V3A 
$} A DT2 -4,-5 cr 
$ } GL0AD cr 
BGLOAD V2A 
$} BGRAPH,BDATA esc 

4, At the ionosonde, turn on the power supply for the line drivers and 
the monitor oscilloscope. Turn the RECORDING START/ INTERVAL SECECTOR 
switch to manual, turn the camera switch off and turn the data 
oscilloscope power supply off. Run one sweep to determine Fof2 and 
then set the frequency as desired. Reduce the transmitter power to 
save tube life as far as is consistant with good echo amplitude. 

5. To start the collection programs, set console switch 0 to 1 and 
switch 17 to 0 and type the following. 

OPERATION 
£ cr 
HEIGHT 

APPROX ECHO HEIGHT cr 

The approximate height is determined from the monitor scope. The 
TTY will print the echoes within 75 km of the height selected. The 
shape of the first echo printed is displayed on the Hazeltine. When 
the exact height has been determined, set console switch 0 to 0 and 
wait until OPERATION is typed on the TTY. Set switch 0 to 1 and type 



the following. 
OPERATION 


1 ^ cr 
HEIGHT 

ACCURATE ECHO HEIGHT cr 
FILE NAME 
XXXXX XXXX cr 

Where the X’s can be any printing characters or spaces such as the date 
and time. The program will type out a heading. Write in the operating fre- 
quency and the critical frequency. To stop the data collection, set switch 
0 to 0. 

6. To transfer the data stored on the disk to paper tape, type the 
following on the foreground TTY. 

+B+C 

FKM15 V3A 

$} A PPO 4/DK2 5 cr 

$}PIP cr 

PIPBF VI 4 A 

> T PPCK-DK2 XXXXXX XXX (A) cr 

To transfer the data to magnetic type, type the following. 
tBtC 

FKM15 V3A 

$} A DT2 4/DK2 5 cr 

$}PIP cr 

PIPBF V14A 

>T DT2+DK2 XXXXXX XXX (A) cr 
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XXXXXX XXX are the same 9 characters typed in for the file name in 
step 5 but with the space after the 6th character instead of the 5th. 



on O O O O O O O O O O O O O O O O O O O O O O O O Q o o 
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PROGRAM LISTING 


C RPROC FIXED FREQUENCY DAT A PROCESSING PROGRAM 

REAL ilDAi'A ( 600 ) , El i < 1 0 ) , ! IK b ) , FN AV E ( 2 ) , S B ( ?. ) 

I NIECE R CDAIA < 5 30 > , f ( 1 02 4 ) , E A < 10 > , I A < 5 ) 

INTEGER TO , T 1 ,T2 , F,F I , PR, OP, 0(5, COL, ALT, DT 


ARRAYS 

CDATA CONTAINS THE AVERAGE OP 600 FRAMES OF DATA 
T IS USED To CONVERT ADC NUMBERS TO Li I NARY 
HD AT A CONTAINS THE HEIGHT Ur THE MAIN ECHO PICKED BY PCK 
FOR EACH 10 SEC INTERVAL 

EH AND EA CONTAIN THE HEIGHT AND AMPLITUDE OF UP TO 10 
ECHOES DETECTED BY ECHO 

Hr AND IA CONTAIN THE HEIGHT AND AMPLITUDE OF UP TO 3 
ONE-HOP AND 2 TWO-HOP ECHOES THAT ARE PRINTED AND 
STORED ON THE DISK 

FNAME IS THE NAME Or THE FILE USED TO STORE HEIGHT AND 
AMPLITUDE DATA. THE FORMAT FUR IMPUTING FRAME IS XXXXX XXXX 
WHERE X IS ANY PRINTING CHARACTER OR SPACE 
SB CONTAINS 11 »• AND "*•» , A * IS PRINTED WHEN MORE ECHOES 
ARE DETECTED THAN CAN BE PRINTED 

VARIABLES 

TO IS THE THRESHOLD USED BY ECHO 
Tl AND T2 ABE THE ONE-HOP AND TWO-HOP THRESHOLDS 
FI IS THE NUMBER OF FRAMES IGNORED DOE TO NOISE OR DUE 
TO ADC TIMING ERRORS 

dG IS THE HEIGHT OFFSET DUE TO THE 30 UScC ENCODE GATE DELAY 
DT IS THE .DAT ASSIGNMENT OF A DISK 
COL IS USED TO PRINT 

ALT IS USED T'J END A LINE WITH AN ALTMODE WHICH SUPPRESSES 

THE NORMAL CARNAGE RETURN 

CO MM ON /D AT A S /C 0 AT A , BO , TO , F R I 

COMMON/DFS/T, HDATA 

DATA SB( I ),Sd(2)/!H , I H*/ 

FS=600 
DT=2 
BG=30 
IFLAG = I 

ALT=( (7*3+6) *3+4) *-51 2 
CQL=( 3*d+5 ) *4096 
FRJ=75000. /I 6334. 

SEf UP TABLE TO CONVERT ADC NUMBERS TO 
INTEGERS IN THE RANGE 21 d GO -2 IB 
DO 1 1-1,512 
X=FLOAT( I-l )*32./75. 

J= I NT ( X ) 

IF (X-FLOAT(J) .GT. .5) J=J+I 


1 T ( I ) =J 

DO 2 1=51 3, 1024 
X=FLUAT( 1 025-1 ) *32. /75. 

J=INT(X ) 

IF (X-FLOAT(J) .GT. .5) J=J+I 

2 T ( I }=-J 

100 WRITE (6, 50) 

50 FORM AT ( 1 OH 1 OPERAT I ON ) 

C INITIALIZE THRESHOLD USED IN SUBROUTINE ECHO 

Tl =10000 
T2 = 1 0000 
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C DATA FDR BACKGROUND PROGRAM 

CALL BGDAK I FLAG, fC:CK,CDAf A, riDATA, Nil, IV,I3,F,T1 ) 
«EAD(6,5J) DP 

C OPERATION = 0 PRINT ALL ECHOES FOUND NEAR XHT 

C = I LUOK FOR ECHOES NEAR XHT AND 2*XHT 


C 

C 

C 

51 

52 

53 


54 

55 

C 

C 

C 

c 

56 


C 

3 

C— 

C 

200 

c 

c 

c 

c 


c 


c 


c 

/> 

c 

c 


c 

c 


PRINT ECHO HEIGHT AND AMPLITUDES AND STORE DATA 
ON .DAT 2 THIS DEVICE SHOULD BE A DISK 
A DECTAPE TAKES TOO LONG 

FORMAT (14) 

WRITE (6, 52) 

FORMAT ( 7H HEIGHT) 

READ(6, 53) XHT 
FORMAT (F5 . 1 ) 

IF COP .EG. 0) Gn TO 3 
WRITE(6, 54 ) 

FORMAT (I OH FILE NAVE) 

READ( 6, 55 ) FRAME ( I ),FNAME<2) 

FORMAT (A5, I X, A4) 

CALL ENTER (DT, FNAME ) 

ONE EXTRA LINE MUST BE ADDED TO THE BEGINNING OF THE' FILE 
AND TWO AT THE END FUR THE HP 9830 GRAPH PROGRAM TO 
WORK PROPERLY 

WRITE(DT,64) 0,600,0,0,0. ,0,0. ,0 
WRITE HEADING 
WRITE (6, 56) 

FORMAT ( 1 X//1 X , 4HDATE»20X , 1 OHFREQ (MHZ) , 20X.9HCRIT FREQ/// 

1 1 5X , 5H I HUP,28X,5H2 HGP/22H TIME FR HGT AMP, 

224X, 9HHGT AMP) 

ARGUEMENTS FOR RDATA 
CALL RDAT(T,CDATA,FI) 

NH=0 


BEGINNING OF DATA COLLECTION 
CALL GO ( LL ) 

LL=I IF SWITCH 0 IS SET 
-1 IF SWITCH 0 IS RESET 
IF (LL .EO. 0) GO TO 20 

NH COUNTS # OF 10 SECOND PERIODS SINCE START 
NH-NH+ I 

IT IS RESET AFTER 100 MIN 
IF (NH .GT. 600) NH=| 

IM AND IS ARE THE TIME IN MINUTES AND 10'S OF SECONDS 

I M = ( NH- 1 )/6 

IS=NH-I-IV*6 

WAIT FOR ADC INTERRUPT ROUT I HE TO FILL CO AT A 
CALL RDATA 
F=FR-F I 

CALCULATE FACTOR TO CONVERT NUMBERS IN CDAT A TJ MIITIVOITS 
FR I “FR J/F LO AT ( F > 

IF (UP .EO. I) GO TO 30 


THE FOLLOWING SECTION PRINTS ECiiU HEIGHTS FOUNT) WEAR 

XHT. IT IS USED TO GET AN ACCURATE HEIGHT FOR STARTING TIT- 

DATA COLLECTION 

IX- 1 NT ( X 1 IT*2 . / 3 . ) - BG+ 1 

TQ«TI 

IX IS THE CENTER OF THE REGION IN CD AT A THAT IS SEARCHED 
FUR ECHOES 

CALL EC1 1’J (IX, Nt: , EH , E A ) 
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II =T0 

IF (NE .EQ. 0) 00 TO M 
C. LOOK FOR THE LARGEST ECHO 

L=l 

DU 10 1=1, NE 

10 IF (EA(I) .GT. EAC L ) ) L«I 

C CENTER THE BACKGROUND ECHO SHAPE GRAPH OH THE HIGHEST ECHO 

KKK = I NT (EH(L)*2 . Z3.+.5) -BG+ 1 
HD AT A ( liH)=EH(L) 

IFLAG=0 

WRITE (6, 57 ) ri f COL,IS,F, ( El!( I ) , I =1 , ME ) 

57 FORMAT! tX, 12, A1 ,11 , 1 HO , 1 4 , 1 0F 6 . 1 ) 

GO TO 200 

C NO ECHOES FOUND 

I | KKK — I X 

C START THE BACKGROUND ECHO SHAPE PROGRAM, DO NOT FIT PARABOLA 

IFLAG=-I 

4 WRITE (6 ,57 ) IM,COL,IS,F 

GO TO 200 

C IF OPERATION = J, CLOSE DATA FILE 

20 IF (OP .EQ. 0) GO TO 100 

I=-l 

WRITE(DT,64) 1,600,0,0,0. ,0,0.,0 
KR I TE ( DT, 64 ) I , 600,0,0,0. , 0 , 0 . , 0 
CALL CLQSE(UT) 

GO TO 100 

C DATA COLLECTION 

30 IX=lNT(XNT*2./3.)-!3G+l 

C WRITE FIRST PART OF LIKE. ALT SUPPRESSES CARNAGE RETURN 

WRITE (6, 58 > I u , CO!., I S , F J ALT 

58 FORMAT ( 1 X , 12 , A1 , 1 1 , 1 HO, 14, I X , A 1 ) 

C ZERO OUT HEIGHT AMD AMPLITUDE ARRAYS 

DO 31 1=1 ,5 
HT< I ) =0 . 

31 I A ( I ) = 0 
HDAfACNH)=0. 

LL=0 

TO=fl 

C SEARCH FOR ECHJES NEAR XHT 

CALL ECHO < I X , NE , EH , EA ) 

Tl -TO 

IF (NE .EQ. 0) GO TO 46 
C FIND ECHO CLOSEST TO XHT 

CALL PCK(J, EH, XHT, ME) 

K = l 

IF (J .EO. 0) GO TO 32 
C MAIN ECHO FOUND 

X=EH( J ) 

HT( 1 )=X 
I Ad ) =HA ( J ) 

EA( J ) =- 1 000 

:<K K = I NT ( X *2 . /3 . + . 5 ) - !M3+ I 
C NE1 XHT 

XHT=X 

C SAVE HEIGHT FOR HEIGHT VS. TIME GRAPH 

HDATA(NH) =X 
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C START BACKGROUND GRAPH PROGRAM 

IFLAG=0 
K-2 

C ARRANGE OTHER ECHOES If! ORDER OF THERE AMPLITUDES 

32 U'J 04 LL=K, 0 

I? (ME .LT. LL) GO TO 35 
1 = 1 

DO 30 J=t , NE 

03 IF (EACJ) .GT. EA(I>) I=J 

i IT ( LL ) =EiK I ) 

I A ( LL)=EA( I ) 

04 £A( I ) =- 1 000 
a =4 

05 LL=LL-1 

C PRINT OUT ONE-HOP ECHOES AND THEIR AMPLITUDES 

C IH+ SUPPRESSES LINE FEED AT THE START OF A LINE 

DO 06 1 = 1 , LL 

36 WRITE (6,59) HTC I ) , I A( I ) , ALT 

59 FOR MAT < l H+ , 1 X , r 6 . 2 , 1 4 , A 1 > 

37 J=2 

C IF THERE ARE MORE THAN 3 ONE-HOP ECHOES, PRINT A * 

IF (ME .LE. 3) J=l 
rip I TE (6,60) S ( J ( J ) , ALT 

60 FORMAT ( 1 M+,2A1 ) 

I X = I NT ( X HT*4 . /3 . ) - EG+ 1 

L=4 

TO =T2 

C LOOK FOR TWO-HOP ECHOES 

CALL ECHlX IX, NE2 , EH, EA) 

T2=T0 

IF (NE2 .EO. 0) GO TO 41 
C SPACE OVER TO fNU-HGP COLUMNS 

L=( 3-LL ) * 1 I 

WRITE (6-, 61 ) ( S3 ( I ) , 1 = 1 ,L) .ALT 

61 FORMAT ( I H+,50AI ) 

C ARRANGE TWO-HOP ECHOES IN ORDER OF THEIR AMPLITUDES 

DO 39 L=4,S 

IF (NE2 .LT. L-3) GO TU 40 
1 = 1 

DU 48 J=I,ME2 

48 IF (EACJ) .GT. EA ( I ) ) I=J 

HT(L)=EH( I) 

IA(L)=EA( I) 

C PRINT uur TWO-HOP ECHO HEIGHT AND AMPLITUDE 

WRITE (6, 59) HT ( L) , I A ( L ) , ALT 

39 EA( I )=-l 000 
L=6 

40 IF (NE2 .LE. 2) GO TO 41 

C END LINE WITH A * IF THERE ARE MORE THAN 2 TWO- HOP ECHOES 

WRITE (6,62) 

62 FORMAT (2H+*> 

GO TO 42 

41 NRITE<6,63) 

63 FORMAT ( 2 H+ ) 

C ELIMINATE SPACES IN ARRAY IH AND HT 

42 L=L-4 

DO 43 1=1 ,L 
J=LL+I 
K=3+ 1 
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riT<j)=;rrm 

43 IA(J)=IA(\) 

L“L+LL * 

IF (L • Cl: • *3) GiJ 10 43 
J=L+1 

00 44 I=J,3 
HT ( J ) =0 . 

44 I A ( J ) =0 

C STORE DATA UN .DAT 2 

45 WRITE (DT, 64) Nil, F , L, LL, C.-il* C I ) , I A< I > , 1= 1 , 2 ) 

64 FORMAT ( !X, 13,14, 12, 12,2 (F7.2, 14 )) 

C IF THERE ARE M'JRE THAN 2 ECHOES TOTAL, STORE THE REST 

C ‘JN THE NEXT LIME 

IF (L ,LE. 2) GO TO 200 

HR ITE( DI, 6b ) ( HT (I), I A (I), I =3,5 ) 

65 FORMAT ( 1 X, 3(F7 .2, 14) ) 

GO TO 200 

C NO ECHOES FOUND, GRAPH AMPLITUDE HEAR XHT B'JT DC NOT 

C FIT PARABOLA 

46 x;<k=ix 

IFLAG=-I 
GO TO 37 
END 


.TITLE RDATA /MACRO SUBPROGRAMS FOR RRPROC 

.GLOSL RDATA, RDAT, 3GDAT, . DA , GU 

.DEC 

NS A? 1=5 30 

l'H=50 

FR=oOO 

.OCT 

CP=204 
SCR=1 12 
TT=6 

GO 0 /SET SH TO I IF DATA SNITCH 0 IS SET 

JMS* .DA / 0 IF DATA SNITCH 0 IS RESET 

JMP .+2 

SH 0 

DZM*’ SH 

LAS 

SPA 

ISZ*- SH 

JMP* GO 

RDAT 0 /PARAMETERS FOR DATA COLLECTION 

JMS* .DA 

JMP .+4 

T 0 /LOG OF LJC OF ARRAY TG CONVERT ADC #'S TO BIN **$ 

BUF3 0 /LOC OF LUC OF AVERAGED INTEGER DATA 

I 0 /LOC OF * OF FRAMES IGNORED 

. INIT TT, I , RDAT 

DZM* (CP /FOREGROUND IS A C CONFIRMER 

LAC* T 

DAC T /PUT LOC OF ARRAY T IN T 

LAC* 3UF3 /PUT LUC OF ARRAY CD AT A IN 8UF 3 
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DAC BUF3 

LAC (NS AM* 3 /SET UP CODS TO TAKE ADC DATA FROM 30F! 

PAL /CUNVEi'iT IT TO INTEGERS IN THE RANGE 21 1 

CLX /TO -218 AND ADD IT TU BUF3 

LAC T2 /STRAIGHT LINE CUDE IS NEEDED TO ACHIEVE HIGH 

A DAC T2 t X /SPEED 

I AC 

AXS 3 

JMP A 

CLX 

LAC T2+1 

B DAC T2+ I , X 

I AC 

AXS 3 

JMP B 

CLX 

LAC T2+2 

C DAC T2+2 ,X 

I AC 

AXS 3 

JMP C 

LAC (NS AM /ZERO OUT BUF2 

PAL 

CT X 

D DZM 3UF2 , X 

AXS I 

JMP D 

LAC (-FR /SETUP NUMBER OF FRAMES TO BE AVERAGED 

DAC F 

LAC (f /SET FLG TO I, SET TO 0 WHEN 600 FRAMES HAVE 

DAC FLG /BEEN AVERAGED 

DZM FIG /ZERO OUT * OF FRAMES DISCARDED AS TOO NOISY 

JMS ADGET /START A TO D CONVERTER 

JMP* RDAT 

BGDAT 0 /DATA FOR BACKGROUND PROGRAMS 

JMS* .DA 

JMP .+12 

FLAG 0 /LOG OF WORD SET TO 0 WHEN AN ECHO HAS BFFN FOUND 

K 0 /CENTER OF MAIN ECHO 

CD AT A 0 /ARRAY OF AVERAGED DATA 

HDATA 0 /ARRAY OF HEIGHT VS TIME DATA 

NH 0 /LOC OF NUMBER OF 10 SECOND PFRIOOS 

IM 0 /LOC OF TIME IN MI NUTS 

IS 0 /LOC OF TIME IN SECONDS 

FRI 0 /LOC OF AMPLITUDE NORMALIZATION FACTOR 

FUR 0 /LG C OF I HOP THRESHOLD 

LAC* CDATA 

DAC CDATA 

LAC* HDATA 

DAC HDATA 

LAC (FLAG /LOC SCR IS UNUSED REG IS TOR IN .SCUM TABLE 

DAC* (SCR /BGD PROG LOOKS THERE TU FIND LOC OF FI AG 

JMP* BGDAT /AND OTHER DATA 

R°ATA 0 /WAIT FOR ALL FRAMES TO BE COLLECTED 

LAC F 

PAX 

LAC FLG 
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TINT 

G 

H 


DIE 

DEAD 


5ZA 

JMP 

.WAIT 

•WRITE 

JMP 

PXA 

TAD 

SMA 

JMP 

DAC 

-T.IMER 
• IDLE 
0 

DZM 

LAC 

SZA 


XT /ALL 600 FRAMES COLLECTED BEFORE PROGRAM GOT TO 

TT,2,MSG2,6 /THIS POINT, SLOW PROC, TIME SCALE IS OFF 

LJ 

/ALL 600 FRAMES HOT COLLECTED YET 
(10 /IF THERE ARE LESS THAN 10 FRAMES LEFT TO Bt 
/AVERAGED IN, GO TO 3 AN!) LOOP UN TILL ALL 
G /HAVE BEEN AVERAGED 

.+4 /OTHERWISE START TIMER FOR (F-1O/60 SECONDS 

0,TINT,0 /AND GIVE TIME TO BACKGROUND 


TINT 

FLO 


/RETURN HERE WHEN TIME ELAPSES 

/SHOULD TAKE 1/6 SEC TO FINISH COLLECTION 
/IF TIME IS UP CONUINUE 


/ 

MS G 1 


M.SG2 


JMP 

J 

/OTHERWISE GO TO J 

ISZ 

FLG 

/RESET FLAG 

DZM 

DIE 

/ZERO CUT DIE, ADC V;i 

LAC 

TFIG 

/TRANSFER FRAMES ION- 

DAC* 

I 


JMP* 

RDATA 

/RETURN 

LAC 

(-1 


DAC 

DIE 


•TIMER 

60, DEAD, 0 /START TIMER 

LAC 

SZA 

FLG 

/LOOP UNTILL ALL FRA 

JMP 

L 


JMP 

0 

H 


0 


/IF 1 SEC TIMER HAS i 

ISZ 

DIE 

/COLLECTED, DIE=0, R 

SKP 


/IF ALL FRAMES HAVE ' 

JMP 

M 

/THE ADC FAILED TO S' 

•RLXIT 

DEAD 


DZM 

DEAD 


LAC 

(NS AM 

/PREPARE TO RESTART 

PAL 

CLX 

/ZERO 

GUT BUF2 

DZM 

BUF2,: 

X 

AXS 

1 


JMP 

N 


LAC 

(-FR 


DAC 

F 


LAC 

(1 


DAC 

FLG 


DZM 

FIG 


JMS 

ADGET 

/RESTART ADC 

• WAIT 

IT 

/PRINT ERROR MESSAGE 

•WRITE 

TT,2, 

‘•'SGI , 12 

JMP 

SW+ 1 

/RETURN TO RPROC 


2000 

0 

•ASCII 

2000 

0 

•ASCII 


/ERROR MESSAGES 
/ ADC 0IED'<15> 

'SP'< I 5> 
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/ ADC INTERFACE INSTRUCTIONS 

ADA' I =703 724 /WRITE INITIAL 1 2 = 

ADSU=70370I /SKIP ON OVERFLOW 

ADST=703/2I /SKIP OH TIMING ERROR 

ADCO*703704 /CLEAR OVERFLOW 

ADCT=703/44 /CLEAR TIMING ERROR FLAG 

• SETUP- 1 35 

• SCUM- 1 00 
AD7jCR=26 
ADCAH-ADNCR+ I 

ADGET d /ADC SERVICE ROUTINE 

JMP ADFST /FIRST TIME THRU GO TO ADFST 

DAC* ( ADA'CR ) /PUT NUMBER OF SAMPLES TO BE COLLECTED IN LOC 21 

LAC (BUFI-1 /PUT LoC OF 3UFI IN LOC 26 

DAC* (ADCAR) 

ADWI /START ADC 

JMP* ADGET 

ADFST LAC* (.SETUP) 

DAC 

JMS* .-I 

ADSO 

AD I NT 

LAC* (.SC0M+5I 

REALTP DAC 

LAC* (.5C0M+D4 

CALL4 DAC 

LAW -NS AM 

DAC ADGET + I 

JMP ADGST+2 

ADINT 0 /ADC INTSRUPT ROUTINE 

/PROCESSING STARTS HERE AFTER ONE FRAME HAS 
/BEEN COLLECTED. THE DATA IN A FR\ME IS 
/ADDED TO DUF2 AND EVERY 600 FRAMES, BUF2 IS 
/TRANSFEREE) TO uUF3 AMD !3UF2 IS ZEROED u'JT 
/THEN Ti IE ADC IS RESTARTED. ALL IMIS MUST EE 
/DUNE IN LESS TI1AM 11.4 MSEC TO PREVENT 
/MISSING A FRAME OR STARTING THE 
/ADC IN THE MIDDLE OF A FRAME 

DBA 

DAC ADS A 

ADST 

SKP 

ISZ ERR /SET ERR TO 1 IF A TIMING ERROR HAS OCCURED 
ADCU /IN THE ADC INTERFACE 

ADC f 

LAC (600000+RTS /SET UP ENTRY TO THE PROCESSING 

JMS* REALTP /SUBROUTINE AT A LONER PRIORITY LEVEL 

JMS* CALL4 

ADINT 
0 

LAC ADS A 

DUR 

JMP* ADINT 

RTS 0 

DAC AREG /SAVE THE A, X AND L REGISTERS 

PX A 

DAC XREG 

FLA 
DAC 


LREG 
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LAC ERR 

UZM ERR /ZERO CUT ERR 

SZA /IF THERE A AS A TIMING ERROR, GO TO S 

JMP S /DATA FOR THIS FRAME IS LOST 

LAC r /SET UP FOR CONVERT I. MG ADC NUMBERS TO INTEGERS 

PAX 

LAC* BUFI f X /IF THE AMPLITUDE OF THE FIRST SAMPLE IS GREATER 

TAD (-TH /THAN A THRESHOLD THE A HOLE FRAME 15 ASSUMED Tu 

SMA /8E NOISY AND IT IS IGNORED 

JMP 3 

12 LAC* SUFI , X /CONVERT NT'/SER IN BUFI TO INTEGER 

TAO S1JF2 /ADD CORRESPONDING ENTRY OF 6UF2 

DAC UUF2 /STORE- SUN BACK IN BUF.2 

.BLOCK 3*N SAM-3 /REPEAT FOR EACH SAMPLE IN SUFI 

P ISZ F /DONE KITH ONE FRAME, ADD ONE TO F 

J’lp R /IF ALL FRAMES ARE nOT COLLECTED, GO TO R 

LAC FIG /ALL FRAMES COLLECTED, SAVE THE NUMBER OF FRAMES 

DAC TFIG /IGNORED 

DZM FIG 

LAC <-FR /SETUP TO RESTART ADC 

DAC r 

DZM FLG 

LAC ( NS AM /TRANSFER BUF2 TO BUF3 AND ZERO OUT BUF2 

PAL 

CLX 

0 LAC BUF2,X 

DAC* BUF3.X 

DZM BUF2,X 

AXS I 

JMP Q 

LAS /IF DATA SWITCH 0 IS SET GO TO R, CONTINUP 

SPA /PROCESSING 

JMP R 

LAC LREG /IF DATA Si'.’ITCH 0 IS 0, LOAD REGISTORS WITH 

PAL /THEIR ORIGINAL CONTENTS AND RETURN 

LAC XREG 

PAX 

LAC AREG 

•RLXIT RTS 

R LAC XREG /RESTORE REGISTERS, RESTART ADC AND RETURN 

PAX 

LAC LREG 

PAL 

JUS ADGET 

LAC AREG 

. RLX IT RTS 

S ISZ FIG 

JMP P 

ADS A 0 

TFIG 0 

F 0 

FLG 0 

FIG 0 

AREG 0 

XREG 0 

LREG 0 

ERR 0 

BUFI .SLOCK NS AM 

BUF2 .BLOCK NS AM 

.END 
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SUBROUTINE ECHG< IX , ME, H;i, EA) 

SUBROUTINE ECHO SEARCHES THE AREA FROM SO DIVISIONS 
BELOW IX TO SO DIVISIONS AMOVE IX FOR ECHOES 
A THRESHOLD IS CALCULATED AS REIMS 5 TIMES THE AVERAGE 
NOISE. AH ECHO IS 4 (JR MORE CONTINUOUS SAMPLES ABOVE 
THE THRESHOLD 
REAL Eli CIO) 

INTEGER CD ( 530 ) , CD AT A ( 1 0 1 ) , £A< I C ) , T, 3G 

COMMON /OATAS/ CD,DG,T,FRI 

NE=0 

M I N =CD ( I X ) 

U=0. 

J=0 

DO NOT SEARCH OUTSIDE OF CD AT A 

IF (IX .LT. SI) IX =5 1 

IF (IX .ST. 475) I X=475 

DO I 1=1,101 

II-IX-5I+I 

K=CD( II) 

AN AMPLITUDE ABOVE THE THRESHOLD IS ASSUMED TO BE 

PART OF AN ECHO AMD NOT INCLUDED IN THE NOISE 

IF ( K .GT. T) GO TO I 

IF (K .LT. MIN) M I N — K 

SUM UP NOISE SAMPLES 

U=U+FLOAT(K) 

J=J+1 

CDATA(I )=K 

IF THRESHOLD IS TOO LOW THERE WILL NOT BE ENOUGH NOISE SAMPLES 
IF (J .LT. 50) GO TO II 

CALCULATE AVERAGE NOISE AND NEW THRESHOLD 
U=U/FLOAT(J) 

T=M I N+5 * ( I NT ( !J ) + 1 -M I N ) 

1 = 1 

SKIP OVER SAMPLES ABOVE THRESHOLD AT BOTTOM OF RANGE 
ALL OF THAT ECHO MAT NOT BE IN THE ARRAY 
IF (CDATACI) .LT. T> GO TO 4 
1 = 1 + 1 

IF (I .GT. 101) RETURN 
GO TO 3 

LOOK FOR AN AMPLITUDE ABOVE THE THRESHOLD 
IF (CDATAU) .GT. T> GO TO 5 
1 = 1 + 1 

IF ( J .GT. 101) RETURN 
GO TO 4 

FIND WIDTH UF REGION ABOVE THRESHOLD AND MAXIMUM AMPLITUDE 
JC— I 

bo 6 J= I ,101 
M=J+ I 

IF ( M .GT. 101) RETURN 

IF (CDATA(M) .GT. CD AT A ( < ) ) K=M 

IF (CDATA(M) .LT. T) GO TO 7 

CONTINUE 

I»M 

IF WIDTH IS LESS THAN 4 DIVISIONS, IT IS NOT AN ECHO 
IF (J CUd 4) GO TO 4 
I = I-J 
NE=NE+ I 

IF (NE .GT. 10) RETURN 


106 


C FIND PARAMETERS OF ECHO AND FIT A PARABOLA TO IT 

Z=CDATA(K-I ) 

C=CDATA(K ) 

C RECORD ECHO AMPLITUDE 

E A < M E ) =C*F E I + . 5 
X-CDATA < K+ 1 ) 

A= ( X-C-C+Z ) /2 . 

B=(X-Z >/2. 

C CALCULATE INTERPOLATED HEIGHT 

Z=-3/<A+A) 

A=A*.6 

C RECORD ECHO HEIGHT 

EH(NE)-(Z+F LOAT ( I X+K+BG-5 2 ) ) *3 . /2 . 

C CALCULATE WIDTH OF ECHO AT A POINT HALF NAY BETWEEN THE 

C MINIMUM AND THE MAXIMUM AMPLITUDE OF THE ECHO 

W= ( ii*B/A- ( 2 . * { C-FLO AT <M I N ) *FR I ) ) ) /A 
IF (W .LE. 0.) GO TO 10 
W=SQRT(W) 

u=K-nrr<w) 

I2=K+INT<W) 

IF (II .LT. I) 11=1 
IF (12 .GT. 101) 12=101 

C SUBTRACT PARABOLA PROM PART OF ECHO ABOVE THE HALF HAY 

C POINT AND SUBTRACT AN INVERTED PARABOLA FROM THE PART 

C OF THE ECHO BELOW THE HALF WAY POINT 

DO B J=I 1 ,12 
X 

Y=(A*X+B)*X+C 

IF (A5S(X> .LT. W/2.) CO TO 3 
U=N 

IF (X .LT. 0.) 0=-N 

Y=- ( A* ( X-0 ) + 3 ) * ( X-0 ) +FLU AT ( M I M ) 

3 CD ATA ( J ) CDATA ( J ) - 1 NT ( V ) + M I N 

IF (NE .EQ. 1 ) GO TO 4 
X=EH(NE)-EH(NE-I ) 

IF (X .LT. 0. ) GO TO 4 

IF (X .GT. 14.) GO TO 4 

X *F LOAT ( E A ( N E ) ) /FLOAT ( E A ( NE- 1 >) 

IF (X .GT. .225) GO TO 4 
10 N£=NE-1 

I =M 

GO TO 4 

C IF THRESHOLD IS TOO LOW, RECOMPUTE IT STARTING WITH 

C A HIGHER VALUE 

J! T-250./FHI 

J=0 
U=0. 

DO 12 1=1,101 

IF ( CDATA ( I ) .GT. T) GO TO 12 
U=U+FLU AT (CDATA ( I ) ) 

J=J+I 

12 CONTINUE 

IF (J .LT. 50) RETURN 

GO TO 13 

END 
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SUBROUTI --4E PCK ( J , H, XHT, NH ) 

SUBROUTINE PCK FINDS THE ECHO CLOSEST TO XHT AND SETS J 
EQUAL TO THE POSITION OF THAT ECHO IN ARR\Y H. IF NO ECHO 
IS WITHIN 20 KM OF XHT, J IS SET TO 0 UN RETURN 
SEAL H (10) 

D=20. 

J=0 

DG 1 1=1 , ME 
A=ABS(tl( I )-XHT) 

IF '(A .CT. D) GO TU I 

J=I 

D=A 

CONTINUE 

RETURN 

END 


BGRAPH 

PROGRAM BGRAPH PRODUCES A GRAPH OF THE HEIGHT VARIATION OF THE 
MAIN ECHO OVER AN II MINUTE PERIOD OR A GRAPH OF THE ONt- 
AND TWO-HOP ECHO SHAPES 
REAL A(25, 3) ,F<o7) 

INTEGER B(25) ,0(25) ,TI < 12) , I3G,FR,T,C0L, Fr 

CO MMOi J/DAT/A , F 

DATA BLAM, STAR/I H , I \\*/ 

C0L=< 3*3+5 )*40D6 

CALL I N IT < B , D , ;C , FR , T , F F ) 

READ SETTING JF SWITCH 17, IF SET, GO TU 2 AND GRAPH HEIGHT 
VS. TIME IF IT IS RESET, TO TO 10 AND GRAPH THE ECHO SHAPE 
CALL SvVR 
I F ( i< ) 1,10,2 

READ STARTING TIME FOR HEIGHT VS. TIME GRAPH 
READ! 6,50) If 
FORMAT ( 12 ) 

CLEAR SCREEN 
CALL Z A P ( 2 S ) 

F‘1AX=0. 

FMI M=l 000. 

FIND PRESENT TIME 
CALL MFRM 
J=IT*o+ I 

IF STARTING TIME IS GREATER THAN PRESENT TIME GO SACK AND 
READ SWITCH SETTING AGAIN 
IF <K .LH. J) GO TO I 

SET END TIME AT K 0!? II HN GREATER THAN STARTING TIME, 
WHICHEVER IS S <ALLER 
IF (K .CT. J+A6) K= J+66 
DO 3 I =J , K 
CALL HCETCC, I ) 

IF (C .ED. 0.) GO TO 3 

FIND MINIMUM AMD ’U.XIMUM HEIGHT UVER THE .RANGE OF THE GRAPH 
IF (C .CT. FMAX ) F"AX=C 
IF (C .LT. FMIN) F'<I.H=C 
CONTINUE 

CALCULATE NORMALIZATION FACTOR 
SP -~ < F M AX— i- M I N ) /2 B . 



OOUl-siOOtfiOtfi O 4* 
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C DIVIDE THE HEIGHT RANGE INTO 25 PARTS 

DU 6 11=1 ,25 
FMIN=F J 'AX-SP 
FH=FMAX-SP/2. 

>t - ] 

DU 4 1=1 ,67 
F(I) = !5LAN 

FIND ECHOES IN PRESENT PART OF HEIGHT RANGE 
D'J 5 I=J, 

CALL HGET(C.I) 

IF (C .EG. 0.) GO TJ 5 
IF (C .Gf. FNAX ) GO TO 5 
IF (C .LE. FMIl!) GO TO 5 
M=I+l -J 
F ( M ) =STAI? 

CONTINUE 

OUTPUT ONE LINE OF THE GRAPH 
WRITE (6, 51 ) Fii, < F ( I ) , I = I , M ) 

FORMAT ( 1 X , F 5 . 1 ,6? A I ) 

FMAX=FMIN 
CONTINUE 

SETUP AND WRITE GUT THE TIME AXIS 
DU 7 1=1,12 
TI ( I ) = IT+I-I 

WRITE ( 6 , 52 ) <T K I ) , 1= 1 , I 2) 

FORMAT (IX, 1216) 

GO TO I 


EC HU SHAPE GRAPH 

1 0 83=30 
FF= 1 

C' FILL ARRAYS B AND D WITH PROPER SECTION OF CD AT A 

CALL 3 DAT A 
N =2 + FF 

C CALCULATE FACTOR TO CONVERT INTEGERS IN B AND D TO MILLIVOLTS 

FRI=75000./< I 6384 . *FLDAT ( FR ) ) 

THR =FLO AT ( T ) *F R I 

C CONVERT DATA IN B AND D TO MILLIVOLTS AND PUT IN REAL ARRAY 

DO II 1=1,25 
J=26- I 

A<J,2)=FLOAT(0(I ))*FRI 
A(J, 1 )=FL0AT(3( I ) >-;.-FRI 

11 CONTINUE 
PM=0. 

C SKIP PARABOLA FIT IF NO MAIM ECHO V.'AS FOUND 

IF (FF , EO. 0) GO TO 12 
CALL PARA (A, P 1 ' ) 

Ki<=K+ 3G-I 

PM=1 . 5 * ( F LO AT ( KK ) + PM ) 

C GET TIME AND CLEAR SCREEN 

12 CALL TI.VE( I", IS) 

CALL ZAP( 28) 

CALL GRAPH (i‘i,KK, I.M,CGL, IS, PM, THR) 

GO TO I 
END 
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.TITLE 

BOAT A 


/ 

THESE a 

;uaRn utiles take data fhum. the foregrgund pitoguamo a* 

/ 

NEEDED 

n dev: 

z LDP T; f LI GRAPHS Uri "TIE BACKGROUND TER?- f i:lAL 


.GLGdL 

SNR , 3D AT A , HGET , I H IT, . DA , NF RM , T I M E 

3CR=i 1 2 




I NIT 

0. 


/INITIAL I Z AT I UN RUUf I NE 


JMS* 

.DA 

/ARGUE! 'ENTS ARE SET UP NON TO SAVE TIME LATER 


JMP 

.+7 


9 

0 



D 

K 

0 

0 



FR 

0 



THH 

0 



FF 

0 




LAC* 

(SCR 

/LOCATION SCR IS AN UNUSED REGISTER IN THE 


DAC 

LUC 

/TABLE. IT CONTAINS THE LOC UF THE ARRAY IN 


LAC* 

LHC 

/WHICH THE LOC OF DATA FOR 8GRAPH IS KEPT 


DAC 

FLAG 

/FLAG IS LOC OF WORD WHICH IS SET TO 0 OR -1 


ISZ 

LOC 

/DEPENDING ON WHETHER AN ECHO WAS FOUND OR NOT 


LAC* 

LUC 



DAC 

;<a 

/KA IS LOC OF WORD CONTAINING INDEX OF CENTER 


ISZ 

LOC 

/OF MAIN ECHO IN CDATA 


LAC* 

LOC 



DAC 

CDATA 

/CDATA CONTAINS LOC OF ARRAY CDATA 


ISZ 

LOC 



LAC* 

LOC 



DAC 

HD AT A 

/HD AT A CONTAINS LOC OF ARRAY HD AT A 


ISZ 

LOC 



LAC* 

LOC 



DAC 

NH 

/LOC OF NUMBER OF 10 SEC. PERIODS SINCE START 


ISZ 

LOC 

/OF DATA COLLECTION 


LAC* 

LOC 



DAC 

LIM 

/LIM AND LIS CONTAIN LOCATION OF TIME IN 


ISZ 

LOC 

/MINUTES AND SECONDS SINCE THE START 


LAC* 

LOC 

/OF DATA COLLECTION 


DAC 

LIS 



ISZ 

LOC 



LAC* 

LOC 



DAC 

F 

/LOC OF WORD CONTAINING OF FRAMES AVERAGED 


ISZ 

LOC 



LAC* 

LOG 



DAC 

T 

/LOC OF WORD CONTAINING ONE-HOP THRESHOLD 


JMP* 

I NIT 



.DEC 



jf 

BDATA 

0 

/THIS 

SUBROUTINE TRANSFERS DATA FROM CDATA TO B AMD D 

A 

LAC* 

FLAG 

/STAY IN THIS LOOP UNTIL FLAG IS SET TO 0 OR - 


SZA1SMA 



JMP 

A 



I AC 




SNA 


/IF FLAG IS -1,. PUT 0 IN FF, NO MAIN ECHO 


DZM* 

FF 



LAC 

(1 

/RESET FLAG TO 1 


DAC* 

FLAG 



LAC 

(-25 

/PREPARE TO TRANSFER 25 WORDS FROM CDATA TO B 


DAC 

COUNT 

/AMPLITUDE NEAR XHT 


LAC* 

B 



DAC 

BB 

/LOC QF ARRAY B 


LAC* 

KA 

/GET CENTER OF ONE-HOP ECHO 
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DAC* 

K 


TAD 

CDATA 


TAD 

(-1 3 


DAC 

r 

C 

LAC* 

i 


DAC* 

SB 


ISZ 

i 


ISZ 

BB 


ISZ 

COUNT 


BMP' 

c 


LAC 

(-25 


DAC 

COUNT 


LAC* 

D 


DAC 

BB 


LAC* 

;<a 


TAD* 

K A 


TAD 

CDATA 


TAD 

(16 


DAC 

I 

CT 

LAC* 

i 


DAC* 

BB 


ISZ 

i 


ISZ 

BB 


ISZ 

COUNT 


JMP 

CT 


LAC* 

1 


DAC* 

THR 


LAC* 

F 


DAC Jr 

FR 

7 

JMP* 

3DATA 

HGET 

0 

JM3* 

.DA 


JMP 

.+3 

CC 

0 


II 

0 



LAC 

HD AT A 


TAD* 

II 


TAD* 

rr 


TAD 

(-2 


DAC 

II 


LAC* 

1 1 


DAC* 

CC 


ISZ 

II 


ISZ 

CC 


LAC* 

II 


DAC* 

CC 


JMP* 

HGET 

/ 


— . 

IIFRM 

0 

LAC* 

NH 


DAC* 

K 

/ 

JMP* 

NFRK 

SWR 

0 



LAS 



AMD 

< 1 


DAC* 

K 


JMP* 

SWR 




/FItfbT WORD OF B IS AMPLITUDE 12 01 VIS IONS BELOW 
/THE CENTER OF THE ONE-HOP ECHO 
/TRANSFER DATA 


/SET UP COUNT FOR SECOND CURVE ON GRAPH. 
/AMPLITUDE NEAR 2*XHT 
/LOC OF ARRAY D 


/TRANSFER DATA 


/GET THRESHOLD AND NUMBER OF FRAMES COLLECTED 


/THIS ROUTINE RETURNS THE HEIGHT OF TH 17 MAIN 
/ECHO FOUND AT TIME IT* 1 0 SECONDS 


/THIS SUBROUTINE RETURNS THE NU’NT-.H OF 10 SECOND 
/PERIODS THAT HAVE ELAPSED SINCE THE START OF 
/THE COLLECTION 


/THIS SUBROUTINE RETURNS I IF SWITCH 17 IS 
/SET AND 0 IF IT IS RESET 



Ill 


TIME 

0 



JMS* 

.DA 


JMP 

.+3 

IM 

0 


IS 

0 



LAC* 

L IM 


DAC* 

IM 


LAO 

LIS 


DAC* 

IS 


JMP* 

TIME 

LTV 

0 


LIS 

0 


CD AT A 

0 


KA 

0 


FLAG 

0 


I 

0 


COUNT 

0 


NH 

0 


HD AT A 

0 


36 

0 


LOC 

0 


r 

0 


F 

0 

.END 



/THIS SUBROUTINE RETURNS TilE TIME IN MINUTES AhU 
/SECONDS SINCE THE S1‘ ANT IF COLLECTION 


SUBROUTINE PARAC A, PM ) 

C SUBROUTINE PARA FITS A PARABOLA TO THE ONE-HOP ECHO AMD 

C DETERMINES THE ASSUMED SHAPE OF THE ECHO 

REAL AC 2-3,3), MIN 
J=l 

C' SET AC 1,3) TO -1000 SO UNUSED PARTS HILL NOT APPEAR 

C UN THE GRAPH 

DO I 1=1 ,25 
AC I t 3)=-i 000. 

I IF CAC I , I ) . LT. AC J, I )) J=I 

C DETERMINE THE PARAMETERS OF THE PARA30LA 

HIN-ACJ, I ) 

C Z IS THE AMPLITUDE 1.5 X-‘ BELO'.f THE CENTER OF THE ECHO 

C X IS THE AMPLITUDE 1.5 KM ABOVE THE CENTER OF THE ECHO 

Z=AC I 4 , 1 ) 

C=AC 1 3, 1 ) 

X=AC 12,1 ) 

AOCZ-'C-C+X )/2. 

B=CX-Z)/2. 

PM=-j/( AC+AC ) 

AC = AC*. 6 

11=3 CRT C C u*F. /AC- ( 2 . * C C-M IN) ) ) /AC ) 

11 = 1 3-IilTCN) 

12=1 3+ I MX (A) 

IF Cl I .LT. 1 ) 11=1. 

IF (12 .ST. 25) 12=23 

C CALCULATE SHAPE OF MODIFIED PARABOLA 

DO 3 1=11,12 
X=1 3-1 

Y=(AC*X + i3)*X+C 

IF CA'SSCX) .LT. . 1 / 2 ,) GO TO 3 
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o=/; 

IF (X .LT. 0.) U=-/i 
Y =- ( AC* ( X-Q ) +0 ) * < X-D ) + M 1 1 i 
8 A ( 1 , 3 ) ~Y 

RETURN 
EMU 


SUBiKMT I N E GRAPH ( .'! , K , I V , COL , IS, PM , T! :R ) 

c THIS SUBROUTINE PRODUCES A GRAPH OF THE AVERAGE AMPLITUDE 

C NEAR XHT AND 2*XHT. IF AM ECHO 15 FOUND NEAR XHT, N = 3 AND A 

C MODIFIED PARABOLA IS FITTED TO THE ECHO AMD DISPLAYED ON THE 

C GRAPH, OTHERWISE H=2 AND MO PARABOLA IS GRAPHED. 

C K IS THE LOCATION OF THE MAIN ONE-HOP ECHO IN CD AT A 

C IM AND IS ARE THE TIME IN MINUTES AND SECONDS 

C COL IS THE CODE FOR THE SYMBOL * 

C PM- IS THE INTERPOLATED HEIGHT OF THE MAIM ECHO 

C THR IS THE THRESHOLD USED BY SUBROUTINE ECHO FOR THE 

C ONE-HOP ECHO 

REAL A(25,3) ,D(4) ,F(67 ) 

INTEGER COL 
CO 'AVON /DAT/ A, F 

DATA D ( I ) , D(2 ) , DC 3 ) , DC 4 )/l HO, I HM , t H+, I H*/ 

DATA DASH, BLAN, EYE, ZERO, f/1 H-, I H , I HI , I HO, t HT/ 

FMAX=-BOO. 

FMI N=500. 

C FIND MINIMUM AND MAXIMUM IN A C I , J ) DO NOT COUNT THE PARABOLA 

C IN FINDING THE MINIMUM 

DO I 1=1 ,25 

IF < AC 1 , 1 > .LT. F M I N ) FMIN=A(I,i) 

IF (AC 1,2) .LT. FMI N) FMIN=A( 1,2) 

DO 1 J=1 , N 

1 IF ( A( I , J) .GT. FMAX ) FMAX=A(I,J) 

IF(FMIN) 2,3,3 

C NEGATIVE MINIMUM 

C CALCULATE 0 POSITION 

2 I Y= I NT ( -F M I N*6 . / ( F MAX- F K I N ) ) + J 
IF (IY .GT. 3 ) GO TO 20 

C CALCULATE AMPLITUDE NORMALIZATION FACTOR 

SP=FMAX/FL0ATC6-IY>*. I 
GO TO. 5 

C LIMIT 0 POSITION OF MIN IS TOO NEGATIVE 

20 I Y=3 

SP— FMIN/30. 

GO TO 5 

C POSITIVE MINIMUM 

3 I Y=0 
SP=F.MAX/60. 

C SET UP FIRST LINE 

5 DO 13 1=2,40 

13 F( I )=UASH 

C PUT AN "1 11 IN EVERY IOTH POSITION 

DO 14 1=1,41,10 

14 F ( I )=EYE 
JZ=IY*I 0+1 

C INDICATE 0 POSITION 

F ( JZ)=ZERQ 
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C WHITE FIRST LINE INCLUDING TIME AND HEIGHT OF ONE-HOP ECHO 

Write (6, 4 ) o-< j),j=i , 4i ), in, col, is, pm 
4 FORMAT! 1 2H ii(Ki) AMP ,41 A) , 14, Al , 1 1 , I HO. 3H I,F3.3,2H 1) 

C DO FOLLOWING FOR EACH OF THE 25 LINES OF THE GRAPH 

DO 16 1=1 ,25 

C MAX INDICATES THE LENGTH OF THE LINE 

C IA IS AMPLITUDE OF ONE-HOP TRACE 

I A= I NT ( A ( I , I I + .3) 

C Fii IS THE CORRESPONDING HEIGHT 

FH-FLOATOC+I3-I )*1 .5 
C BLANK OUT THE LINE 

DO 7 J= I ,61 

7 F(J)=i3LAN 

C ADD 0 AMPLITUDE LINE 

F( JZ)=EYS 

C ADD SYMBOLS FOR EACH OF THE OTHER TRACES 

DO 3 J=I,N 

IM=INT(A( I, J>/SP+.5)+JZ 
IF (IM ,LT. 0) GO TU 3 
IF ( I 'I .GT. MAX) HAX = IM 
IF(FUM)-BLAN) 9,10,9 

C IF THERE IS ALREADY A SYMBOL IN THAT PLACE, PUT A * THERE 

9 F ( IM) =D( 4 ) 

GO TO 8 

10 F ( I M ) =D ( J ) 

8 CONTINUE 

6 IF (I .NE. I) GO TO 16 

C INDICATE THRESHOLD ON THE THIRD LIME 


J=IMT (THR/5P+ .5 ) +JZ 
IF (J .LT. I) J=61 
IF (J .GT.61) J =6 1 
IF (J .GT. MAX ) MAX=J 
F ( J ) =T 

16 WRITE(6 , II) FH , I A, (F ( J ) , J= 1 , MAX ) 

1 I FORMAT! I X , F 5 . 1 , 15, IX, 61 Al ) 

RETURN 

END 



.TITLE 

ZAP 


.clobl 

ZAP, .DA 


. IODEV 

6 

ZAP 

0 



JM5* 

.DA 


JMP 

.+2 

ARG 

0 



.WAIT 

6 


LAN 

(-1000 


DAC 

CNT 

LOP 

ISZ 

CNT 


JMP 

LOP 


LAC 

<176. 


JUS 

PIN IT 


LAN 

(-1000 


DAC 

CUT 

BACK 

ISZ 

CNT 
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JMP yACK 

LAC* ARC 

JMS PR hi 

JMP* ZAP 

PRNT 0 

DAC CUT 

.WRITE 6,3, BUF, 3 
.3 A IT 6 

J.MP* PRNT 

EUF 2000 

0 

JUT 0 

0 

CNT 0 

.END 
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APPENDIX II. OPERATING INSTRUCTIONS AND A LISTING OF PROGRAM GRAPH 

The following procedure is used to graph the data from the digital 

ionosonde once the data is on paper tape. A listing of the graph program 

is included after the operating instructions. 

1. Turn on the calculator and the plotter. Insert the tape cassette 
marked Ionosonde Programs, type L0AD1 and press EXECUTE. 

2. Adjust the plotter to the size of paper used and put a black pen in 
the plotter. 

3. Turn on the paper tape reader, insert the paper tape and press READ. 

4. At the calculator, press RUN and EXECUTE. Enter the. number of 10 sec 

periods of data on the tape. If 0 is entered, the entire tape is 

read. The data are stored in the magnetic tape casette in blocks of 
20 min. 

5. When the data are on the magnetic tape, type in a heading for the graph 
such as the date and time of the data. If several graphs of the same 
data are needed, after one has been made, type CONT 520 and press 
EXECUTE to return to this point. 

6. Enter the limits of the portion of the data to be graphed in minutes 
and seconds. If 0, 0 is typed in as the START and END times, all the 
data are graphed. 

7. Type in the minimum height to be graphed. All echoes below this height 

are ignored. The calculator now makes one pass through the tape to 

determine the minimum and maximum heights and then plots an axis and 
the heading. 

When the display reads PEN? change the pen to a blue one, type 0 and 


8 . 
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press EXECUTE. The height data are now plotted for one one-hop and 
one two-hop echo. 

9. When the display again reads PEN? change the pen to a red one, type 0 and 
press EXECUTE. The one-hop and two-hop amplitudes are now plotted. 

10. When the display again reads PEN? put the blue pen back in the plotter, 
press 0 and EXECUTE. The extra one-hop and two-hop echoes are now 
added to the graph . 
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10 COM USE 5 j 1. 20 1 ? FI I E 5 ? 120] 
20 DIM O I E 4 ] j V$[ 30 ] 


5@ 

60 

70 

n r o n * > r 12 * 

0 

MO " * 


[J J. i ! k 

INPUT 16 

Vlo It 


80 

IF 16 >0 

TUCk! -1 

1 1 J L., 1 'J .!■ 

00 

90 




100 

IF Nl+I 

7*600 

>= 16 THEN 460 

1 1 0 

ENTER C 

1 5 120) 

NlsFlsEsEls (F0RI=lTQ2sHC I P M II :> BE 

1 2& 

FORMAT 

F 3 „ 0 it F 

4 „ 0 it 2 F 2 .i 8 ? F 7 . 2 $ F 4 « @ it F 7 « 2 ? F 4 ,« U 

1. .”! f : l 

IF N1--0 

THEN 

100 

140 

IF N 1 < 0 

THEN 

460 

150 

IF N1=N 

THEN 

220 

160 

IF BBS ( 

N+4-N1 

) < 4 THEN 398 

i 7 8 

I 7= I 7+ 1 



1.30 

N=N“600 



1 90 

IF N1=N 

THEN 

9 c\ 

■i<ii Ii(n ‘nr 

200 

IF ABS< 

N + 4-N 1 

) < 4 THEN 390 

2 1 0 

STOP 


► 

y y \}\ 

IF E < 3 

THEN 2 

50 

2 3 0 

ENTER < 

1 s 240 } 

(FOR I =3105:. HE I » M]> AC I s M It) 

O ft 

FORMAT 

F7 . 2 > F 

4. 0:. F 7 „ 2s F4. 0s F7. 2sF4.0 

■>«. i.. 1 

' » ! I-'I 

p y ..... p p* | 

ir n d t ■*: 

! <J 1 \ .1. 

TO 1 

STEP -1 


280 IF E2< 1-3 THEM 330 


3 8 0 

ML i |I !’i 

3=HC J , M 3 

3 1 0 

Fit: 3: » m ]=ac J s n it 

320 

GOT 0 

360 

338 

hi: i: 5 m 

] “ A E I s M ]=0 

340 

goto 

360 

330 

IF I> 

El THEN 33 

360 

NEXT 

I 

378 

G08UI3 

1 230 

2 8 Pi 

GOTO 

3. 0 0 

29 Pi 

FOR I 

~1 TO 5 

400 

HE 1 1 M 

+ 1 ]--HE I s M 1 

4 1 8 

AC IsM+1 ]=AC Is Ml 

420 

HE I s M 

3= AC I s M 3=0 

430 

NEXT 

I 

448 

GO SUB 

1 230 

450 

GOTO 

150 

460 

FOR I 

= 1 TO 5 

478 

FOR J 

=M TO 120 

480 

HE I s J 

3 = AE I s J 3=- 

490 

NEXT 

;.j 

500 

NEXT 

I 

510 

STORE 

DATA F 

520 

F I HI) 


530 

I) ISP 

" NAME " 5 


M ]) 
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540 INPUT Y $ 

550 rasp "START " 5 
560 INPUT 11,12 
570 1 1 = 1 1*6+12/10+1 
580 laSP "END"! 

590 INPUT 12,13 
600 12=12*6+13^10+1 

610 rasp "min"; 

620 INPUT 1 19 

638 IF 1 2 > 1 THEN 658 

640 12=1439 

650 15=1 

660 GOSUB 1340 

670 Ol 1 ]=0E 3 >0 

680 DC 2 3= Of: 4 '!= 1 000 
690 IF HC1,M]<0 THEN 340 
7 U 0 P U h! J = 1 I U 2 
710 K=J + J 

720 FOR I=3#J-2 TO K+l 
730 S=HCI,M] 

740 IF S=0 THEN 790 
750 IF SCOCK-l] THEN 770 
760 OC K-l ]=S 

770 IF S >OC K 1 OR S<M9 THEN 790 
780 OC K ]=S 
790 NEXT I 
800 NEXT J 

810 IF N=I2 THEN 850 
820 GOSUB 1230 
830 GOTO 690 

840 I 2=N 1 

850 FIND 3 

860 SCALE -0 .1,1. 05 , - 1 „ 25, 1 „ 25 
870 FOR J=1 TO 2 
880 K= : J + J 


890 8=0 1 K 

i :hj[ i< i 

900 OFFSET 

0, 1 . 2*J-2. 3 

910 GOSUB 

;[ 3 3 0 

920 NEXT J 


930 LABEL 

K 'it o Ji 2 !• 0 ;i 0 (< h ) 

940 PLOT 0 

.. 5, 1 „ 05, 1 

950 CPLOT 

-LEN(YT)/2,0 

960 LABEL 

( * ) Y $ 

970 LABEL 

( * , 1 . 7 , 2 ? 0 , 0 . 

9 WO PUR Q= 

1 TO 2 

990 GOSUB 

2070 

1000 FOR L 

= 1 TO 2 

1010 R = 3*| 

O 

1020 OFFSE 

T 0, 1 , 2*L-2. 3 

1030 PLOT 

1 i; 1 ”7 1 X| 

0 . 85 1 j. . 07 , 1 

.1. U ,- r f \ L. 1 1... 

1050 IF Q= 

2 THEN 1110 

1060 LABEL 

<*> "HEIGHT" ! 
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1070 s= O [ K-l !3 — OC K ] 

1080 I3=R 

1.090 GOSUB 1698 

3.100 GOTO 1150 

1110 CP LOT 0J--6.87 

1120 L Fl B E: L ( * > " Fi M P L I T U D E ” ! 

1130 13-R+l 

1140 GOSUB 1698 

1150 NEXT L 

1160 NEXT Q 

1170 GOSUB 2070 

1180 OFFSET 0J-1.1 

1190 1 3=3 

1200 GOSUB 1698 

1210 PLOT 1.05)2.35) 1 

1220 STOP 

1230 N =4-1+1 

1240 M=M+1 

1250 IF N< 12 1 THEN 1380 
1260 M=1 

1270 IF 15=1 THEN 1310 
1280 STORE D Fl T Fl F 
1290 F'”F+ 1 
1300 RETURN 
1310 PEN 

1320 LORD DATA F 

1330 GOTO 1290 

1348 F== I NT < < 1 1-1 >/120::'+3 

1350 N-I 1 

1360 M=I 1 “■ < F-3 > * 1 20 
1378 GOTO 1318 
1380 I3=<I2-Il>/6 
1398 LABEL (*? 1 . 7? 2? 0? 0. 6> 
1400 FOR 1=0 TO 13 STEP 10 
1410 PLOT I / 13 j 0)1 
1428 CP LOT -3.3»-l.l 
1430 IF Kite THEN 1450 

1 4 4 0 C P L 0 T 0 ,, 4 * I N T ( 1 + 1 / 1 0 0 ) j 0 
1458 LABEL < 1470) 1 5 

1460 NEXT I 
1470 FORMAT F4.0 
1480 X=20 

1498 IF S/X>4 THEN 1510 
1580 N=10 

1510 FOR 1=0 TO S/X 
1520 PLOT 6:> I+X/S) 1 

1 5 3 0 C f : ' L 0 T - 6 . 5 1 - 0 . 3 

1 5 4 9 L A B E L. < 1 5 6 0 > ( 0 C K 3 + 1 * X > 5 

1550 NEXT I 

1560 FORMAT F5» 0s " . " 

1570 XAXI8 0? 10/13)0? 1 
1580 VAX IS 0 ) / 8 ) 0 ) 1 
1590 PLOT -0.055)0.3) 1 
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1 6 0 0 L A B E L < * » 1 . 7 - 2 , p I / 2 r 0 . 6 > " H E I G H I < K M ) " 5 
1610 LABEL 1 ,, 7 > 2 ? 0-0.6) 

1.620 PLOT 0 . 5 j 0 > 1 

1630 CPLCi'l" -6 f 2.1 

1640 LABEL <*> “TIME <MIN>"3 
1650 PLOT 0, 2.- 1 „ 02? 1 
1668 LABEL C16?0)J 
1670 FORMAT F3.0.-" HOP" 

1680 RETURN 

1690 GO SUB 1340 

1700 Y"0 

1710 P=I2-I1 

1720 IF 13=3 THEN 1800 

1730 J=2 

1740 X=HC Rt M 3 

1750 IF X< OL K :i THEN i860 

1760 IF ABS<X-Y)>3 THEN 1780 

1770 GOTO 1880 

1780 PEN 

1 8 0 0 G O 8 U B 1 3 0 F 3. 9 0 0 - 1 9 2 0 > 1 9 4 8 j 1 9 0 0 ? 1 9 2 0 
1810 PLOT CN- 1 1 >/Pj Gj J 

1 8 2 0 I i x i 

1830 IF N™ 12 THEN 1880 

1840 GOSUG 1230 

1858 GOTO 1720 

I860 J=1 

1870 GOTO 1810 

1880 PEN 

1890 RETURN 

1900 G=<X-0[ l< K'/S 

1910 RETURN 

1920 G=AL R !> I'l 1/1000 

1938 RETURN 

1940 FOR J ::: 1 TO 2 

1950 K==J + J 

I960 S=OC K-l 3-OC K 3 

1970 FOR I =K+J-1 TO kf+1 

1988 X=HC I ii M 3 

1990 IF X <= M9 THEN 2830 

2 0 0 0 G == ; < X - 0 1 K 3 ) / 8 + 1 . 2 * J - 1 , 2 
2010 PLOT (N-ID/PjG 

2020 PEN 

2830 NEXT I 

2840 NEXT J 

2050 J=1 

2960 RETURN 

2070 PLOT 1 , 05 ? 1 . 15)1 

2 O 8 O F 3. N l.i 3 

2090 D1SP 

2100 DISP "PEN" 3 

2110 INPUT 13 

2 1 20 RETURN 
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APPENDIX III. OPERATING INSTRUCTIONS AND A LISTING OF PROGRAM SHAPE 

Program SHAPE calculates the virtual height and amplitude for a 
reflector and provides easy methods for modifying the reflector shape to 
make the calculated virtual height and amplitude match the experimental 
data. 

1. Turn on the calculator and plotter, insert the tape cassette marked 
Ionosonde Programs and type LOADKEYO and press EXECUTE. Put the 
function key overlay marked SHAPE over the function keys . 

2. Press RUN and INIT (f^) . The calculator asks for the limits of the 
reflector size. If a portion of a graph produced by GRAPH is to 

be analyzed, a sheet of graph paper is put on the graph with the X 
axis lined up with the lowest virtual height and the Y axis lined up 
with the first sample point. The X coordinate of the last sample 
point and the Y coordinate of the highest virtual height, both in 
units of grid intervals are entered into the calculator. If the two- 
hop height in grid units is to graphed, the height variation of the 
two-hop echo is entered next. Otherwise go to step 3. 

3. Press POINTS (f^) and enter the coordinates, in grid units, of key 
points such as relative minima and maxima which have the same height 
as the reflector height. 

4. Points are connected with a half cycle of a sine curve by pressing 
HALF-SINE (f^) • Enter the X coordinates of the two points to be con- 
nected with a half sine and press EXECUTE. 

5. To connect two points with a straight line, press CURVE FIT (f^) . 
Enter 2 as the number of points and then enter the X coordinates of 
the two points to be connected with a straight line. A curve through 
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more than two points is made by specifying the number of points to be 
greater than two and then enter that number of X coordinates of points 
to be fitted. The first coordinate is the starting point of the curve 
and the last point entered is the end point of the curve. 

6. FILL Cf-^g) is used to modify the reflector curve shape by adding or 
subtracting a sine-shaped segment from the reflector. The start and 
end points of the segment to be modified and the amount in grid units 
to be added are entered. 

7. Place the graph paper in the plotter and press MOVE PEN (f . Press 
the space bar to move the pen and adjust the pen position controls 
until the pen moves between the origin of the graph and the coordinate 
made up of the X ordinate of the last sample and the Y ordinate of the 
maximum height. The resulting plots will be the same size as the 
graphed data. 

8. Press PARAM (f^) . Enter the estimated horizontal length of the reflec- 
tor in km, the height of the lowest echo and the difference between the 
maximum and minimum virtual height . 

9. Press AXIS (f,) to draw and label an axis for the one-hop echo. Shift 

fg or f^ draws and labels an axis for the two-hop echo. 

10. Press REF SHAPE (f^) to draw the one-hop reflector shape. Shift fj 

or f^y draws the two-hop reflector shape which is twice as high as 
the one-hop reflector at all points. 

11. Press RES HEIGHT (fg) to draw the one-hop virtual height calculated 
from the one-hop reflector. Press f . „ to draw the two-hop virtual 

lo 

height . 

12. Press AMP (fg) to draw the amplitude of the one-hop echo. Press f^g 
to draw the two-hop echo amplitude. 
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13. In the case of multipath, NORM (f^) changes the scale factor used in 
graphing the reflector shape and the virtual height so that the maxi- 
mum excursions in height fit on the graph. It should not be used when 
trying to match a reflector to a data graph because the scale size may 
be changes. Shift f,. or f ^ does the same thing for the two-hop graph. 

14. HEIGHT (fjj^) displays the height of the reflector in grid units for 


any X coordinate. 
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10 

REM K! 
REM I! 

T I T LfiL I 

2 AT I ON 

3@ 

DIM H* 

5C256 3» 

USC7j73jCSC7 

40 

SCREE 

~0 , 08 t 

1 8 6 1 ! 0 ;> 1 !' 3 

50 

LABEL 

( \\ 1 |i 2. 

!l fc) |l 1 .) 

60 

70 

DISP 

“UPPER 

LIMITS Xj Y"5 

80 
9 0 

INPUT 

Y9=Y5 

)<5 V 5 



3» FSC7 3 


100 DISP "VERT SIZE OF 2 HOP Y"5 
110 INPUT Y9 
120 END 


16 REM KEY 1 

20 REM INPUT POINTS 

30 DISP "INPUT POINTS X <■ Y " 5 

40 INPUT X 1 > Y 1 

50 H[ FNR < 0 >3 : =Y 1 /Y5 

60 GOTO 30 


10 REM KEY 2 

20 REM HALF SINE 

30 DISP "HALF SINE FROM XI TO 

40 INPUT X 1 i X2 

50 Y1“H[ FNFK0) ] 

60 Y2=H[ X2 ! 

70 T=P I / ( X2--X 1 > 

80 T 1 - ( Y2+ Y 1 > *0 „ 5 

90 T2=T1-Y2 

100 FOR I=X1 TO X2 

110 HE I ] ~ T 1 + T 2 * C 0 S ( ( I X 1 ) * T > 

120 NEXT I 

130 GOTO 38 


10 REM KEY 3 
20 REM CURVE FIT 
30 DISP "tt OF POINTS"! 

40 INPUT N 

50 REDIM I.JCN« Nil nCt N3» FUN] 

60 FOR 1=1 TO N 

70 DISP "POINT" 31 3 

80 INPUT XI 

90 FII I 3= HI" FNFl < 0 ) 3 

100 UCI>13-1 

110 FOR J=2 TO N 

120 UC I » J ]=UE 1 1 J-l 3*X 1 

130 NEXT J 

140 NEXT I 

150 11=10:. 2 3 

160 I2=UI:Nh21 



0 IIS- < HE I 3 + HC I l-l : 

0 K 1 =RTN < M) 

0 K1=K1+K1 
0 RETURN Z 

I-HC I + 1 

0 REM KEY 11 
0 REM DISPLAY HEIGHTS 
0 HI SR " X 1 " 5 
0 INPUT XI 

0 DISP H[ FNR<e> ]#Y5> 

0 GOTO 3@ 


8 REM KEY 12 
0 REM MOVE PEN TO 0,0 

AND 3. , 

0 PLOT 0 :i 0,1 
O UNIT 32000 
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58 PLOT 1.1,1 
68 WRIT 32600 
70 GOTO 30 


10 REM KEY 13 

20 REM M CODIFY CURVE 

30 n I SP " START j STOP > I NCR I MENT " 5 

40 INPUT X1,X2,D 

50 Z=FNfl ( 0 > 

60 J=X2-X1 

70 D-D/Y5 

80 FOR I- XI TO X2 

9 0 H [ I ] = H [ I ] + D * S I N ( C I - X j. > * p I / J > 

100 NEXT I 

110 GOTO 38 


10 R 

EM KEY 

5 

28 R 

EM NORMALIZE 

30 V 

1 = H0 


48 V 

2-H0+FI 


50 P 

1=0 


Q PI J- 

70 F 

Ihi« | 

OR 1 = 2 

TO 255 

80 2 

-FNB ( 0 ) 


98 V 

-l<*SQR < 

1 + M * M > 

100 

X-- l<*M + (. 

I - 1 ) *P5 

1 1 0 

IF V > V 1 

THEN 130 

1 20 

V 1 = V 


130 

1 d I/I 

IF V< V2 

'i. 1 V ::::: 

THEN 150 

J. u r U 

150 

IF X>P1 

THEN 178 

160 

P1=X 


179 

IF X<P2 

THEN 190 

130 

P 2 = X 


190 

NEXT I 


200 

V3=V2-V 

1 

2 1 0 

P3-R2 f : ' 

1 

220 

END 



10 REM KEY 6 
20 REM AXIS 
30 PLOT 1,0 
40 PLOT 0,0,2 
50 PLOT 8, Ii- 


IO DEF FHFK2) 

20 X1-INT<X 1*255/ 

40 RETURN XI 


1 
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6 @ D I S P " L fl B E L S T E P S I 2 E I N K M X > Y " J 

70 INPUT X 1 f Y 1 

a 0 T ™ I N T ( V 2. / Y 1 )' * V 1 

9 0 P L O I 0 , 0 1 j < I V 1 > / v 3 ? - 2 

100 PLOT 0? ( T-V 1 >/V3? - 1 

110 CP LOT -4.3? 0.3 

1 2 O L fl B E L ( •S' ) T 5 
130 T-T-’Y1 

140 IF T >= VI THEN 90 
150 T=INT<P1/X1>*X1 
160 IF T=P 1 THEN 180 
3.70 T ™T+N 1 

1 8 0 P L 0 T ( T - F : ' 1 ) / P 3 ? 0 0 1 j 2 
190 PLOT (T-Pl )/P3? 0? -1 

200 CP LOT 1 „ 3 ? -- 1 , 3 

210 IF T :s 0 THEN 230 
2 2 0 C P L 0 T - 1 N T < L G T ( T > ) / 2 ? 8 
230 LABEL < * > T !i 
248 T-T+X 1 

258 IF T <« P2 THEN 180 
260 END 


10 REM KEY 7 

20 REM ASSUMED REFLECTOR SHAPE 
38 J — 0 

48 FOR 1 = 1 TO 256 

50 IF CI-1>*X5/2550>J THEN 100 

6 0 P L 0 T < ( I - 1 > * P 5 •- P 1 > / P 3 ? < H i: I ] * fl + 1-1 0 - V 1 } / V r : i « - ? 

70 NEXT I 

80 PEN 

90 END 

100 J = J+1, 

110 PEN 
120 GOTO 60 


10 REM KEY 8 

20 REM RESULTING HEIGHT 

30 J=0 


40 

FOR 1 = 


TO 

255 


50 

Z-FN8 ( 

0 ') 




60 

V=K*SQ 

R< 

1 +M 

*M ) 


70 

X=K*M+ 

<1 

-1) 

*P5 


80 

IF <1- 

1) 

#X5 

•-"2550 >. 

J THEN 130 

90 

PLOT < 

X- 

P 1 ) 

/R3 ? <:v- 

-V 1 ) -"'V3 ? -2 

100 

NEXT 

I 




110 

PEN 





1 20 

END 





130 

JssJ+1 





140 

PEN 





150 

GOTO 

90 
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10 REM KEY 9 

2 0 R E M R E S U L T I N G E C H O f : l M P L I T U D E 
30 Rl = R--< ;[ 0 1'- y 9 > 

40 . J — 0 

50 FOR 1=2 TO 255 
60 Z-FNB (9 ) 

70 V=K*SQR < 1+M*M> 

80 IF M2=0 THEN 230 
90 R= ( < l+M*M)tl . 5). •-'M2 
100 IF R=V THEN 150 

1 1 0 1 F R 1 > 6 Fi N D R 1 < V R N D R > V T H E N 1 5 0 

120 IF R 1 >V AND R < V RNI) R>8 THEN 150 

1 3 0 Fi 3 = H 0 * S 0 R < fi B 8 < R / < R y ) ) ) / y 
140 IF R3< 10 THEM 160 ’ 

150 R3=10 

160 N=K*M-K 1-1 >*P5 

170 I F ( I - 1 > # X 5 •-•' 2 5 5 © > J T H E N 2 5 0 

1 8 0 P L 0 T < X - P 1 > / P 3 j ft 3 •■•" 1 S s 2 

190 R 1 =R 
200 NEXT I 
210 PEN 
228 END 
230 ft 3 “H0. 

240 GOTO 160 
250 J=J+1 
268 PEN 
270 GOTO 180 


10 REM KEY 15 
20 REM NORMALIZE 2 HOP 
30 V4-H0+H0 
40 V5=V4+ft+ft 


50 P 
60 P 
70 F 

4 ■■■ 0 
OR 1 = 

£ i U .j.j 

80 2 

: ~FNB < 

0 ) 

90 I 

F M > 1 

THEN 230 

1 00 

V = l<* ( 

l+l/COS<K 1 > > 

110 

X =K*7 

ftN < K 1 ) + < I 1 ') * 

120 

IF V) 

V4 THEN 140 

138 

V4=V 


148 

IF V< 

V5 THEN 160 

1 58 

n' K [ :™: \ t l 


160 

IF X> 

P4 THEN 180 

170 

P4==X 


1 80 

IF X< 

P8 THEN 280 

190 

f i' , ‘i 


200 

NEXT 

I 

2 1 0 

V6 = V5 

-V4 


220 P7=P8~F'4 
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10 REM KEY 16 
£0 REM 2 HOP AXIS 
30 PLOT 1 j 1 . £5 
40 PLOT 0)1,25)2 
50 PLOT 0 !i 2 ii £ 5 [i 1 

68 DISP " LABEL STEP SIZE IN KM X) Y"> 
70 INPUT' X 1 ) Y 1 
80 T-INT<V5/Yn*Yl 

9 0 P L 0 T 0 ,, 0 1 ) < T -• V 4 > / V 6 * Y 9 / Y 5 + 1 . 2 5 j - 2 

1 0 0 P L 0 T 0 ) ( T -- V 4 > / V 6 * Y 9 / Y 5 + 1 . 2 5 t ■■■ 1 
110 CPLOT -4 . 3 ? -0 , 3 

120 LABEL OO'r; 

130 T-T-Y 1 

140 IF T >= V4 THEN 98 
150 T=:ENTCP4/N1>*N1 
168 IF T-P4 THEN 180 
170 T-T+Xl 

180 PLOT < T P4>/P7) 1 .26) -2 

1 9 0 P L 0 T < T - P 4 > / P 7 ) 1 . 2 5 > - 1 
200 CPLOT -1.3) -1.3 
218 IF T=0 THEN 230 
2 2 0 C P L 0 T - 1 N I < L G T < T ) > / 2 ) 0 
230 LABEL <*)TS 
248 T=T+N 1 

250 IF T <= P8 THEN 180 
260 END 


10 REM KEY 17 

20 REM ASSUMED REFLECTOR SHAPE £ HOP 
38 J 0 

40 FOR 1=1 TO 256 

50 I F < I - 1 > * X 5 / 2 5 5 0 > J T H E N 1 0 8 

6 0 P L 0 T < < I - 1 > * f : ' 5 - p 4 ) / P 7 :■ < 2 * < H C I I * Fl + H 8 > -- V 4 ■' / V 6 * Y 9 

78 NEXT I 

80 PEN 

98 END 

108 J=J+1 

118 PEN 

120 GOTO 60 


5+ 1 25 ) 


10 REM KEY 18 

28 REM RESULTING 2 HOP ECHO PATTERN 
38 J=0 

40 FOR 1=2 TO 255 
50 Z=FNB(0) 

68 IF M > 1 THEN 140 
78 V=K*(l + l/COS<K:l)> 

86 X»K*TAN<Kl> + <I-n*P5 

98 IF <I-1)#X5/2550>J THEN 160 

1 8 8 P L 0 T < X - P 4 ) / P 7 ) < V V 4 > / V 6 * Y 9 / Y 5 + 1 . ? 5 , - ? 

118 NEXT I 
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120 PEN 
130 END 
140 PEN 
150 GOTO 1:16 
1 66 J=J+1 
170 PEN 
180 GOTO 100 


1 M 

RFM KFY 

t 9 

20 

REM RE Sill 

..TING 2 

30 

R=~< 10+9 1 

;i ) 

4 U 

K2 = H0 


50 

W I™ 14= -140 


60 

FOR 1=2 ' 

TO 255 

,w 'i 

r K) 

,M“NB(0) 


80 

S=K/C03 C 1 

< 1 > 

90 

V=K + S 



fill PL. ITU DE 


100 IF 112=0 THEN 270 
110 R=C<l+M*Mm.5)/l12 
120 IF M >1 THEN 290 
130 W-R*$/<S+S~R) 

1 4 0 I = R * < K + K - M > / < K + K + K + K - N - H - R > 

1 50 IF ( ( W > K AND W 1 < l<2 ) OR ( IKK AND Wl>l<2>> THEN : 
160 IF C ( N >K+K AND W1<K2+K2> OR <U<K+K FIND N 1 >K2-H 
1 7 0 fi 3 = H 0 * 8 Q R < R B S < W * T / ( < K + K - U > * ( T S ) ) } ) / S 
180 IF A3< 10 THEN 208 
190 A3= 10 

2 0 0 X = K * T A N ( K 1 > + ( I - 1 ) * P 5 

2 1 8 P L 0 T < X - P 4 > / P 7 > A 3 / 1 8 + 1.1, ~ 2 

228 K2=K 

230 N1~W 

248 NEXT I 

250 PEN 

26B END 

8 7 0 H 3 H 0 'd 'S '6 LI R ( 0 , !;j -Sr 3 / ( \! ) ) 

230 GOTO 200 
290 PEN 
300 GOTO 240 
31© JaJ+1 
320 PEN 


1 90 

C2> ) THEN 190 
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APPENDIX IV. OPERATING INSTRUCTIONS AND A LISTING OF 
THE PROGRAM TO COMPUTE LATERAL DEVIATION 

This program computes the separation of the ordinary and extraordinary 
reflection points using the formula for lateral deviation given in Budden 
[1961]. The separation of the reflection points is used to find the 
direction of propagation of the disturbance as discussed in Section 2.6. 

To run the program, turn on the calculator and printer, insert the 
cassette marked ionosonde programs and type L0AD2 and press execute. 

When the program is loaded press run and execute. Enter the ionosonde 
operating frequency and the critical frequency when the calculator asks 
for them and a table similar to Table 2.1 is printed. 



10 DIM DC 2 ]» HI 2] 

20 FIXED 4 
30 HC 1 ]=HI £ ]=X=0 

40 DISP "OPERATING FREQUENCY IN MHZ ” 5 
50 INPUT F 

60 PRINT 'OPERATING FREQUENCY = "JF 
70 DISP "CRITICAL FREQUENCY IN MHZ" ? 

30 INPUT F2 

90 PRINT "CRITICAL. FREQUENCY ="5F2 

100 F=F*10t6 

110 F2=F2*10tb 

120 I I 1=F2*F 30. 6 

130 21=200 

140 X1=80.6/<F*F> 

150 B=0. 56*10t(-4> 

160 L=-0 , 15-''0. 56 
170 N2=0. 54/0. 56 
1 80 Y = ( 2 . S * 1 0 1 1 6 > * B / F 

190 PRINT " LATERAL DEVIATION DX/DZ Y="!Y!" 1 

260 PRINT " ZCKM) X 0 X ' 0 H 

210 Y1=N2*Y 

220 Y2=L*Y 

230 K=0.05 

240 Z=J=0 

250 FOR L=1 TO 28 

260 Z=Z+K 

270 N = N 1 * < i ~ I. C Z-Zl) / Z 1 > T 2 > 

280 X3=X 

290 X=N*X1 

300 FOR 1=1 TO 2 

310 IF X >= l-(2-I)*Y THEN 420 

320 0=1+1 -3 

330 1 = Y 2 s Y 2 / 2 + ( 1 — X ) ) 

3 40 U = 1 - X / ( i - T +0 * S Q R < T * T + Y 1 * Y 1 > ) 

3 5 0 D I I J = -■ U + Y 1 + Y 2 + (. 1 — ! ! ) / 6 0 R i. 4 + Y 1 + V ] t if 1 — X j ( j — y \ + y 4 4 '• 

360 HE I 3— HC I 3+ DC I 3*K 
376 NEXT I 
38@ NEXT L 

390 WRITE <15>400)Z>HC 1 Is HE 2 3 » DC 1 DDE 2 1 > Ns X 
408 FORMAT F7.2 j4F9.4jE10.2iF9 . 4 
410 GOTO 259 
420 IF J=1 THEN 490 

430 WRITE < 1 5 1 400 > Z j HC 1 3 j HE 2 3 ? IT 1 .1 j DC 2 1 1 N 1 X3 
446 ,1=1 
456 S1=HC 1 1 
460 S2=Z 
470 HC i ]=DE 1 3=6 
430 GOTO 370 
490 IF 1=1 THEN 370 

50 0 WR I T E (15-40 0 > Z 1 H [ 1 3 - H E 2 3 1 D [ 1 3 - D [ 2 3 , H , X 3 
510 X=HE 2 3-S 1 

530 PRINT" 

540 WRITE (15? 550.) Xj Y 

550 FORMAT "LATERAL SEPARATION IS" 1 F7. 3- "KM 
566 STOP 


VERTICAL SEPARATION 



